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It is shown that the a— uj model which is widely used in the study of nuclear relativistic many-body 
problem can exactly be treated as an Abelian massive gauge field theory. The quantization of this 
theory can perfectly be performed by means of the general methods described in the quantum gauge 
field theory. Especially, the local U(l) gauge symmetry of the theory leads to a series of Ward- 
Takahashi identities satisfied by Green's functions and proper vertices. These identities form an 
uniquely correct basis for the renormalization of the theory. The renormalization is carried out in the 
mass-dependent momentum space subtraction scheme and by the renormalization group approach. 



With the aid of the renormalization boundary conditions, the solutions to the renormalization 
group equations are given in definite expressions without any ambiguity and renormalized S-matrix 
elememts are exactly formulated in forms as given in a series of tree diagrams provided that the 
physical parameters are replaced by the running ones. As an illustration of the renormalization 
procedure, the one-loop renormalization is concretely carried out and the results are given in rigorous 
forms which are suitable in the whole energy region. The effect of the one-loop renormalization is 
examined by the two-nucleon elastic scattering. 
Cn ! PACS: ll.10.Gh, 11.15.-q, 13.75.Cs 
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The quantum hadrodynamics (QHD), as a relativistic quantum field theory for baryons and mesons, has been 
widely applied to studying various nuclear phenomena including the hadron-hadron interaction, the hadron-nucleus 



. i scattering, the bulk and single-particle properties of nuclei, etc. [1-5]. It is commonly recognized that although 
the quantum chromo dynamics is a fundamental theory for strong interaction, the QHD, as an effective field theory 
formulated in terms of hadronic degrees of freedom, provides a simple and reliable approach to produce the nuclear 

■ observables that are insensitive to the short-range dynamics. There are various QHD models, renormalizable and 
, "^T* ', nonrenormalizable, which were tested in the past to reproduce the empirical nuclear properties and the experimental 

data. Among these models, the a — uj model proposed by Walecka [1] has been raising particular interest. This 
| £3 ' model contains proton, neutron and isoscalar, Lorentz scalar and vector mesons a and uj and in the tree diagram and 

■ £3 nonrelativistic approximations leads to a nucleon-nucleon interaction potential which behaves as short-range repulsion 

and medium-range attraction. The early development of this model is based on the relativistic mean-field and Hartree 
, • i-h , approximation and shows that the model is quite successful in the applications to the infinite nuclear matter and 
atomic nuclei. Since the model is renormalizable, it is necessary to consider higher order perturbative corrections to 
the results given in the mean field approximation by a certain renormalization procedure. Along this line, a number 
of efforts were made previously [6-19]. Especially, the efforts were mostly concentrated on the renormalization of 
the model in the study of the nuclear matter at finite temperature and density. In this renormalization, the loop 
expansion and spectral function methods were applied to evaluate the loop corrections. However, there are various 
difficulties to occur in the renormalization [8-16]. For example, in Ref. [12], the authors calculated the nuclear matter 
energy density up to the two loop level and found enormous contributions arising from the loop terms that alter 
the description of the nuclear bound state qualitatively. Therefore, it was concluded that "the loop expansion does 
not provide a reliable approximation scheme in renormalizable QHD" [19]. To this end, one may ask what is the 
correct procedure of performing the renormalization for a model of QHD? and how to assess the applicability of a 
renormalizable model of QHD for which the renormalization is carried out? To answer these questions, it is meaningful 
to examine the renormalization of a QHD model from different angles and, as suggested in Ref. [19], "to develop 
and apply systematic and consistent "power counting" schemes that lead to more general conserving approximations 
and to study renormalization group methods that could determine the analytic structure of the ground-state energy 
functional" . 

In this paper, we confine ourself to discussing the renormalization of the a — uj model by the renormalization group 
method in the case of zero temperature. The procedure is very similar to that described in our previous work on 
the QED and QCD renormalizations [20] . The main features of the renormalization given in this paper contain two 
aspects: (1) the renormalization is based on the U(l) gauge symmetry because the a — uj model, as argued in the 
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next section, is exactly of the U(l) local gauge symmetry; (2) The renormalization is carried out by a mass-dependent 
momentum space subtraction [21-24] which will lead to rigorous renormalized results by the renormalization group 
method [25-28]. Ordinarily, the massive vector fields such as the oj meson field, the p meson field and so on are not 
viewed as gauge fields because the mass term in the Lagrangian is not gauge-invariant [29-31]. On the contrary, it was 
pointed out in Refs. [32-34] that a massive vector field must be viewed as a constrained system in the whole space 
of the vector potential (x) . This is because a massive vector meson has only three polarization states which need 
only three spatial components of the vector potential A fl (x) to describe them. While, the remaining component of the 
A fl (x) appears to be a redundant degree of freedom which must be eliminated by introducing the Lorentz condition. 
According to the general principle for constrained systems, the gauge-invariance of a massive Abelian or non-Abelian 
gauge field should be seen from its action given in the physical space defined by the Lorentz condition. This viewpoint 
will be explained in more detail in the next section. From this viewpoint, it is easy to see that the a — ui model is 
surely of U(l) local gauge symmetry Therefore, the model may be quantized by the method as used in the gauge field 
theory In this paper, we will describe the Lorentz-covariant quantization performed in the both of Hamiltonian and 
Lagrangian path-integral formalisms by following the procedure proposed in Refs. [32-34]. From this quantization, we 
obtain an effective action which contains a gauge-fixing term and a ghost term in it and manifests itself to be invariant 
under a set of BRST transformations [35] . It should be mentioned that the quantum theory of the a — u> model was set 
up previously by the method of canonical quantization and in the path-integral formalism [29-31, 36, 37 ]. Especially, 
with the time paths being generalized to a manifestly covariant form, a covariant path-integral formulation for the 
model at finite temperature was achieved in Ref. [37] and led to manifestly covariant Feynman rules for both real 
and imaginary times. Nevertheless, owing to lack of the gauge- fixing term and ghost term in the effective action, the 
generating functional given in these quantizations would not exhibit the BRST-symmetry. 

As emphasized in Ref. [20], a correct renormalization procedure for a gauge field theory must respects the gauge- 
symmetry (the Ward-Takahashi identities [38, 39]), the Lorentz-invariance (the energy-momentum conservation) and 
the mathematical convergence principles. Otherwise, the renormalization would be incorrect. From the gauge- 
invariance (or say, the BRST-symmetry) of the a — uj model, we derive a set of Ward-Takahashi (W-T) identities 
satisfied by the generating functionals, Green's functions and vertices which provide a firm basis for the renormalization 
of the model. As mentioned before, in this paper, the renormalization of the a — uj model will be performed in the 
mass-dependent momentum space subtraction. The prominent advantage of such a subtraction is that it naturally 
provides boundary conditions satisfied by the renormalied wave functions, propagators and proper vertices for the 
quantum a—uj model. These boundary conditions enable us to uniquely determine the solutions to the renormalization 
group equations for those renormalized quantities. With the solutions of the renormalization group equations, a S- 
matrix element can be expressed in the form as given in the tree diagrams provided that the physical parameters in 
the S'-matrix element are replaced by the effective (running) ones. To specify the procedure of the renormalization 
group method, the one-loop effective physical parameters are concretely calculated and given exact and analytical 
expressions. 

The remainder of this paper is arranged as follows. In Sec. II, we present arguments for the gauge-invariance of 
the a — ui model. In Sec. Ill, the a — uj model will be respectively quantized in the Hamiltonian and Lagrangian 
path-integral formalisms. In Sec. IV, we will derive a set of W-T identities obeyed by the generating functionals. In 
Sec. V, a W-T identity satisfied by the uj meson propagator will be derived and the renormalization of the propagator 
will be discussed. In Sec. VI, we will derive a W-T identity satisfied by the vectorial vertex (nucleon-nucleon-w meson 
vertex) and discuss the renormalizations of the vertex and the nucleon propagator. In Sec. VII, the renormalizations 
of the a meson propagator and the scalar coupling vertex (nucleon-nucleon-cr meson vertex) will be derived and 
discussed .Sec. VIII is used to sketch the renormalization group method and the renormalized S-matrix elements. 
Sec. IX serves to derive the one-loop effective coupling constants and masses. In the last section, summary and 
discussions will be made. In Appendix A, the gauge-independence of the S-matrix elements given in the one-loop 
level will be proved. In Appendix B, we will show the differential cross section of the two- nucleon elastic scattering 
in the approximation of order g 2 and examine the effect of the one- loop renormalization on it.. 



II. ARGUMENT OF GAUGE-INVARIANCE FOR THE a - to MODEL 



The a — uj model is described by the following Lagrangian density [1] 




(2.1) 
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the nucleon isospin doublet in which ipp and ip n are the proton and neutron field functions respectively. 



0^ = 0^- ig v A^ - -g s l^ (2.3) 

is the covariant derivative in which and ip stand for the u> and a meson fields, g v and g s designate the vectorial 
and scalar coupling constants, 

= dfnA,, - d v A„ (2.4) 

is the vector field strength and M, m w and m a are the masses of nucleon, u> meson and a meson respectively. In the 
above Lagrangian, the scalar self-couplings are ignored as was done originally in the Walecka model [1]. 

In the previous, the a — u model was considered to be gauge-non-invariant with respect to the following local U(l) 
gauge transformations [29-31] 

A' li (x)=A ll (x)+d li 9(x), ( '> 
<p'(x) = <p(x) 

where 6{x) is the scalar parametric function of U(l) group since the mass term of the u> meson in the Lagrangian is 
not gauge-invariant. But, this does not mean that the dynamics of the lu meson system is not gauge- invariant. As 
mentioned in the Introduction, the oj meson field must be viewed as a constrained system in the space spanned by the 
four-dimensional vector potential A fl (x). As we know, a massive gauge field has three polarization states which need 
only three spatial components of the four-dimensional vector potential A^ to describe them. In the Lorentz-covariant 
formulation, a full vector potential A fl (x) can be split into two Lorentz-covariant parts: the transverse vector potential 
Aj,(x) and the longitudinal vector potential A^(x), 

A^(x)=A^(x)+A l " L (x) (2.6) 

where 

A£(x) = [<r - ^d»)A v {x), (2.7) 



A*{x) = ^Wa,(z) (2.8) 

with □ = d^d^ being the D'Alcmbcrtian operator. The vector potentials A^(x) and A^(x) satisfy the following 
transverse and longitudinal field conditions (identities): 

O^x) = 0, (2.9) 



(g^-^d v )Al(x)=0 (2.10) 

and the orthogonality relation 

J d 4 xA^(x)A Lfl (x) = (2.11) 

which characterizes the linear independence of the two field variables. Since the Lorentz-covariant transverse vector 
potential At^(x) contains three-independent spatial components, it is sufficient to represent the polarization states of 
a massive vector boson. Whereas, the Lorentz-covariant longitudinal vector potential A^ appears to be a redundant 
unphysical variable which must be constrained by introducing the Lorentz condition: 

X = 9"i M = (2.12) 
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whose solution is 



A% = 0. (2.13) 
With this solution, the a — w model Lagrangian may be expressed in terms of the independent dynamical variables 

C = ^(id„ + g v A Til + \g.w) - M]ip - \F^F Tltv + \mlA» T A Tll 

where F^ u is defined as in Eq. (2.4) with replacing the A^fa) by At^(x). The Lagrangian represented above gives a 
complete description of the dynamics of the a — u) model. If we want to represent the dynamics in the whole space 
of the full vector potential as described by the Lagrangian in Eq. (2.1), the w field must be treated as a constrained 
system. In this case, according to the general procedure for constrained systems as formulated in Mechanics, the 
Lorentz condition in Eq. (2.12), as a constraint, must be introduced from the onset and imposed on the Lagrangian 
in Eq. (2.1) so as to guarantee the redundant degree of freedom to be eliminated from the Lagrangian. Otherwise, 
the Lagrangian in Eq. (2.1) itself can not give a complete description for the u> field system. From the Lagrangian in 
Eq. (2.14), one may derive an equation of motion satisfied by the u> meson field as follows: 

d ft F^ + mlA^ = -f (2.15) 

where 



f = 5 ^7^ (2.16) 

is the current generated from the nucleon field. The above equation describes the evolution of the independent variable 
Aj, with time. In particular, when we take divergence of the both sides of Eq. (2.15), considering the identities in 
Eq. (2.9) and d v d^F^ v = 0, we immediately obtain the current conservation 

0*% = (2.17) 

which shows that the current is transverse. 

Ordinarily, the Lorentz condition is viewed as a consequence of the following uo field equation of motion which is 
derived from the Lagrangian in Eq. (2.1) [29-30, 36] 

duF^+mlA" = -j v (2.18) 

The argument of this viewpoint is as follows. When we take divergence of the equation (2.18) and notice the current 
conservation, it is found that 

m*0M„ = (2.19) 

Since m w ^ 0, the above equation leads to the Lorentz condition which implies that one component of the vector 
potential is not independent. It is pointed out here that the above viewpoint actually is an ill-concept and the procedure 
leading to the Lorentz condition logically is not consistent with the principle established well in the mechanics for 
constrained systems. In fact, the aforementioned derivation seems to imply that the Lorentz condition has already 
been included in the Lagrangian denoted in Eq. (2.1). If so, when the Lagrangian is written in the first order form, 
we should see a term in the Lagrangian which is given by incorporating the Lorentz condition with the aid of the 
Lagrange multiplier method. Nevertheless, as will be shown in the next section, there is no such a term to appear 
in the Lagrangian. Moreover, as we know, equations of motion should describe the evolution of the independent 
variables with time as the equation given in Eq. (2.15) does and should not lead to a constraint condition which 
implies some variable in the equation is not independent. Therefore, the viewpoint stated above is not reasonable. In 
accordance with the general principle for constrained systems, the correct procedure is to treat the Lorentz condition 
as a primary constraint and to impose this condition on the Lagrangian in Eq. (2.1) from the beginning. The necessity 
of introducing the Lorentz condition can also be seen from the derivation mentioned in Eqs. (2.18) and (2.19). The 
equation (2.19) can be understood in such a way that if the Lorentz condition is not introduced, there would appear 
a contradiction that the right hand side of the equation is zero, but the left hand side is not. Only when the Lorentz 
condition is introduced, the contradiction disappears. In this case, due to the Lorentz condition, the equation (2.19), 
as a trivial identity, naturally holds and the equation of motion (2.18) can naturally go over to the equation (2.15), 
exhibiting the self-consistency of the theory. Particularly, in the latter case, when the divergence of the equation (2.18) 
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is taken and the Lorentz condition is employed, one immediately obtains the current conservation in Eq. (2.17). In 
addition, we would like to note that for the quantum theory, in the zero-mass limit: m u — > 0, the vector field part of 
the Lagrangian in Eq. (2.1) naturally goes over to the one for the massless vector meson, but, as shown in Sec. V, the 
vector meson propagator does not and occurs a worse singularity, revealing a severe inconsistence of the theory. Only 
when the Lorentz condition is introduced initially and incorporated into the Lagrangian by the Lagrange multiplier 
method, a consistent quantum theory can be constructed. 

Now, let us turn to address the gauge-invariance of the a — u model. Usually, the gauge-invariance is required 
to the Lagrangian. From the dynamical viewpoint, as pointed out in Refs. [32-34], the action is of more essential 
significance than the Lagrangian. This is why in Mechanics and Field Theory, to investigate the dynamical and 
symmetric properties of a system, one always starts from the action of the system. Similarly, when we examine the 
gauge-symmetric property of a field system, in more general, we should also see whether the action for the system is 
gauge-invariant or not. In particular, for a constrained system such as the massive vector field, we should see whether 
or not the action represented in terms of the independent dynamical variables is gauge-invariant. This point of view 
is easy to understand from the mechanics for constrained systems. Suppose a mechanical system is described by a 
Hamiltonian 

H(jp i ,q i )(i = l,2,---,n) (2.20) 

which is given in the 2n-dimensional phase space and constraint conditions 

f a {p ll Qi) = 0(a = 1, 2, • • • , 2m < 2n) (2.21) 

which define a physical phase space of dimension 2(n — m) where the system exists and moves only. If the constrained 
variables can be solved out from the constraint conditions, we may write a Hamiltonian 

H*(p*,q*)(j = 1,2,- -^n-m) (2.22) 

which is expressed via the independent variables and gives a complete formulation of the constrained system. Obvi- 
ously, to examine some symmetry of the constrained system, it is only necessary to see if the Hamiltonian H*(p*,q*) 
other than the Hamiltonian H(pi, q{) to have the desired symmetry because in contrast to the H*(p*,q*), the H(pi, qt) 
is not complete for describing the system. 

Certainly, in some special cases, the Lagrangian given in the physical space itself is locally gauge-invariant so that 
the gauge-invariance of the corresponding action is ensured. This situation happens for the massless gauge fields 
and the massive Abelian gauge field. The gauge transformation of an Abclian gauge field was shown in the third 
equality in Eq. (2.5). Since d li 6{x) acts as a longitudinal field, according to the decomposition denoted in Eq. (2.6) 
and considering the independence of the fields Aj,(x) and A^{x), the gauge transformation of the ui field can be 
equivalently divided into two transformations: 

A%{x) = A^{x) (2.23) 



A'1(x) = A*(x) + d»6{x) (2.24) 

Eqs. (2.23) and (2.24) clearly express the fact that the gauge transformation only changes the unphysical longitudinal 
part of the vector potential, while, the physical transverse vector potential is a gauge-invariant quantity. Furthermore, 
it is easy to verify that the longitudinal vector potential A^{x), which may be expressed as A^ L (x) — d^(p(x) where 
ip(x) is a scalar function, is cancelled in the field strength tensor so that 

F" u = d^A" - d u A" = d»A^ - d" 'A£ = F£ v (2.25) 

This indicates that the longitudinal part of the vector potential has no kinetic energy term in the Lagrangian and 
hence has no any dynamical meaning. Such a vector potential can only be viewed as a constrained variable. Since 
the transverse field variable Aj, is gauge-invariant, the Lagrangian (2.14) which is written in the physical space is 
manifestly gauge- invariant. Therefore, the action given by this Lagrangian is gauge- invariant. Alternatively, the 
gauge-invariance may also be seen from the action given by the Lagrangian in Eq. (2.1) which is now constrained by 
the Lorentz condition. Under the gauge transformation written in Eq. (2.5) and the Lorentz condition denoted in 
Eq. (2.12), it is easy to find that 

SS = -m* j d 4 xed^A„ = (2.26) 
This indicates that the a — v model can surely be set up on the basis of gauge-invariance principle. 
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III. PATH-INTEGRAL QUANTIZATION OF THE o - ui MODEL 



A. Quantization in the Hamiltonian path-integral formalism 

According to the general procedure of dealing with constrained systems, the Lorentz condition (2.12) may be incor- 
porated into the Lagrangian (2.1) by the Lagrange undetermined multiplier method to give a generalized Lagrangian 
[32, 40]. In the first order formalism [32, 40, 41] , this Lagrangian can be written as 

C x = ^(i7"D M - M)V + hd^if - \m\^ + \F^F^ 

where A^ and F M „ are now treated as the mutually independent variables and A is chosen to represent the Lagrange 
multiplier. Using the canonically conjugate variables defined by 

= -j = iV7°, (3-2) 
dtp 

%=^ = 0, (3.3) 

d V 

dC 

n, = - = ^ (3.4) 

and 

n /vi dL j? i \x r-Ffco = Ek, if = fc = 1, 2, 3; . , 

" { ) = 9> = M0 + ^ = { A = -So, if M = 0, (3 ' 5) 

the Lagrangian in Eq. (3.1) may be rewritten in the canonical form 

C = E»A„ + n^i> + Il tp tp + A C-EoX-H (3.6) 
where = (Eq, Ek) is a, Lorentz vector, 

C = d^E^ + m 2 A + gJh°*l>, (3-7) 
X was defined in (2.12) and 7i is the Hamiltonian density expressed by 



U = \{E k f + \{F tl f + \ml [(Aol + [A k f] + jplj + (Vrf 
+m 2 r ip 2 ] - i-ip^f ■ Vip + Mipip - gv^l kr ^Ak - g s ^Tp<p 



(3.8) 



in which i*^ was defined in Eq. (2.4). In the above, the four-dimensional and the spatial indices are respectively 
denoted by the Greek and Latin letters. Eq. (3.6) clearly shows that the terms AqC and Eox are respectively given 
by incorporating the constraint condition 

C = (3.9) 

and the Lorentz condition into the Lagrangian by the Lagrange multiplier method and the Lagrange multipliers Ao 
and E are just the constrained variables themselves in this case. Since the A and E are a pair of the canonically 
conjugate unphysical variables, their constraint conditions in Eqs. (2.12) and (3.9) should simultaneously occur in 
the Lagrangian (3.6). Otherwise, if the Lorentz condition is not introduced, the term E x does not appear in the 
Lagrangian shown in Eq. (3.6) or in Eq. (2.1). In this case, the Lagrangian could not be complete for describing the 
constrained system under consideration. 

From the stationary condition of the action constructed by the Lagrangian (3.6), one may derive the following 
first-order canonical equations of motion: 

A k = d k A Q - E k , (3.10) 
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E k = d l F lk + mlA k + d k E Q + g v ipj k i>, (3.11) 

P = U V , (3.12) 

n v = V 2 p - ml(p + g s i/jip, (3.13) 

{i-fd li -M + g v 'fA li +g a <p)1> = 0, (3.14) 

^(iy^/i + M — grfA,, - g s p) = (3.15) 

as well as the constraint equations written in Eqs. (2.12) and (3.9). Eqs (3.10) and (3.11) act as the equations of 
motion satisfied by the independent canonical variables A k and E k (k = 1,2,3) which precisely describe the three 
degrees of freedom of polarization for the massive u> field, while, Eqs. (2.12) and (3.9) can only be regarded as 
the constraint equations obeyed by the constrained variables Aq and Eq because in these equations, there are no 
time-derivatives of the dynamical variables A k and E k . It is clear to see that in Eqs. (3.10)-(3.15), (2.12) and (3.9), 
there are altogether twelve equations. They are sufficient to determine the twelve variables including the dynamical 
canonical variables tp, vp, W Vl ip, A k and E k (k = 1,2.3) and one pair constrained variables Ao and Eo, showing the 
completeness of the equations. 

Now, we turn to formulate the quantization performed in the Hamiltonian path-integral formalism for the a — u> 
model. In accordance with the general procedure of the quantization, we should first write a generating functional of 
Green's functions in terms of the independent canonical variables which are ip, ip, II V , ip and the transverse parts of 
the vectors A^ and for the u meson field [32, 40, 41] 

Z[J", J,rj, V ] = ±fD(A$.,E£,il>, n v ,<p)exp{i J d 4 x[E»A T(L + . 

+U v ip-n*(A^,E^,tP,tP,U ip ,p) + J£A Tfl + jtp + rM + i/m]} 1 • ' 

where r H*(Aj,,Ej,,tl>,tp,H v ,ip) is the Hamiltonian which is obtained from the Hamiltonian (3.8) by replacing the 
constrained variables A^ and E£ with the solutions of the equations (2.12) and (3.9) 

H*(A»,E£,---)=H(A»,E»,---) \ x=0 ,c=o (3.17) 

and J[t,J,r) and rj are the external sources coupled to the u> meson, a meson and nucleon fields respectively. As 
mentioned before, Eq. (2.12) leads to A^ L = 0. Noticing this solution and the decomposition 

E»(x) =E£(x)+E£(x), (3.18) 

when setting 

E£{x) = d£Q{x) (3.19) 
where Q(x) is a scalar function, one may get from Eq. (3.9) an equation obeyed by the scalar function Q(x) 

a x Q(x) = W(x) (3.20) 

where 

W(x) = -^(zfrVfc) - mlA° T (x). (3.21) 
With the aid of the Green's function G(x — y) (the ghost particle propagator) which satisfies the equation: 

a x G(x-y) = 6 4 (x-y), (3.22) 
one may find the solution to the equation (3.20) as follows: 

Q(x) = ( d 4 yG(x-y)W(y). (3.23) 
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From the expressions given in Eqs. (3.19), (3.21) and (3.23), we see that the E£(x) is a complicated functional of the 
variables Aj, and Ej, so that the Hamiltonian TL* (Aj, , Ej, , • • •) is of a much more complicated functional structure 
which is not convenient for constructing the diagram technique in perturbation theory. Therefore, it is better to 
express the generating functional in Eq. (3.16) in terms of the variables A^ and E^. For this purpose, it is necessary 
to insert the following delta-functional into Eq. (3.16) [32, 40, 41] 

S[Al]S[El - E£(A° T ,^)\ = detMS[C]S[ X ] (3.24) 

where M is the matrix whose elements arc 



M( X ,y) = {<?(*), *(»)} " 
= U x 5\x-y) 



SA^z) 5E»(z) 6E„(z)6A»(z)) (3.25) 



where {C(x), x(u)} i s the Poisson bracket as defined in the second equality in Eq. (3.25). The relation in Eq. (3.24) 
is easily derived from Eq. (2.12) and (3.9) by applying the property of delta-functional. Upon inserting Eq. (3.24) 
into Eq. (3.16) and utilizing the Fourier representation of the delta-functional 

8[C] = J D(p/2^)e L i d ' xp{x)c{x \ (3.26) 

we have 

Z[J", J,%v] = jfjD(A fl ,E^^,U^,ip,p/2Tr)detM5[ X }exp{iJd 4 x^A^ 
+n^ + n v ip + pC-H{A„,E^,ip,i;,n v ,ip) + J»A„ + Jip + vi) + 'il;Ti}}. 1 ■ ' 

In the above exponent, there is a E -ve\ateA term E (doA — dop) which permits us to perform the integration over 
Eq, giving a delta- functional 

S[d A - d oP ] = det|^|-M[A) - P]- (3.28) 

The determinant cZet |c? 1 — 1 , as a constant, may be put in the normalization constant N and the delta- functional 
S[A — p] will disappear when the integration over p is carried out. The integrals over Ek 11^, are of Gaussian- type 
and hence easily calculated. After these computations and noticing the expression in Eq. (3.2), we arrive at 

Z[J»,J,% V ] = ±J D(A^^,ip,)detM_6[d"A^exp{i J d A x[C ( . 

where C was written in Eq. (2.1). When employing the familiar expression [41, 42] 

detM = J D ^^ c yfd*xd*yC( X )M( x ,y)C(y) = J D ^^ c yf ^ 3Q) 

where C(x) and C(x) are the mutually conjugate ghost field variables and the following limit for the Fresnel functional 

J[0MJ - lim C[a}e-^I dix(d " A ^ 2 (3.31) 
where C[a] <~ IlC^m;) 1 ^ 2 an< ^ supplementing the external source terms for the ghost fields, the generating functional 

X 

in equation (3.29) is finally given in the form 

Z[J», J,r?,7?,£,£] = ^jD(A ll7 ip^,(p,C,C)exp{iJd 4 x[C ef f , , 

+JM„ + J<p + rft> + il>ri + ZC + CZ]} 1 ' ' 

where 

jC ef f = C -^( d ^ A ^) 2 + C UC ( 3 - 33 ) 

which is the effective Lagrangian for the quantized a — ui model in which the last two terms are the so-called gauge- 
fixing term and the ghost term, respectively. In Eq. (3.32), the limit a — > is implied. Certainly, the theory may be 
given in arbitrary gauges (a^O). In this case, as will be seen shortly, the ghost particle will acquire a spurious mass 
v = y/am^. 
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B. Quantization in the Lagrangian path-integral formalism 



Now let us to quantize the a — uj model in the (second order) Lagrangian path-integral formalism following the 
procedure proposed in Ref. [32-34, 40]. For later convenience, the Lagrangian in Eq. (2.1) and the Lorentz constraint 
condition in Eq. (2.12) are respectively generalized to the following forms: 

Cx = C - |A 2 (3.34) 

and 

d^A^ + aX = (3.35) 

where A(a;) is an extra function which will be identified with the Lagrange multiplier and a is an arbitrary constant 
playing the role of gauge parameter. According to the general procedure for constrained systems, Eq. (3.35) may be 
incorporated into Eq. (3.34) by the Lagrange multiplier method to give a generalized Lagrangian 

C x = C + Xd^A^ + ^aX 2 . (3.36) 

This Lagrangian is obviously not gauge- invariant. However, for building up a correct gauge field theory, it is necessary 
to require the dynamics of the gauge field to be gauge- invariant. In other words, the action given by the Lagrangian 
(3.36) is required to be invariant under the gauge transformations shown in Eq. (2.5). By this requirement and 
applying the constraint condition (3.35), we have 

5S X = -\J d A xd v A v {x){U x + v 2 )6(x) = (3.37) 

where v 2 — am 2 ,. From equation (3.35) we see ' A v = — A ^ 0. Therefore, to ensure the action to be gauge- 
invariant, the following constraint condition on the gauge group is necessary to be required 

{n x + v 2 )6{x) = 0. (3.38) 

The constraint condition in Eq. (3.38) may also be incorporated into the Lagrangian in Eq. (3.36) by the Lagrange 
undetermined multiplier method. In doing this, it is convenient, as usually done, to introduce the ghost field variable 
C{x) in such a fashion 

6{x) = qC{x) (3.39) 

where <; is an infinitesimal Grassmann's number. Based on the above definition, the constraint condition (3.38) can 
be rewritten as 

(n x + v 2 )C = (3.40) 

where the number <j has been dropped. This constraint condition usually is called ghost equation. When the condition 
(3.40) is incorporated into the Lagrangian (3.36) by the Lagrange multiplier method, we obtain a more generalized 
Lagrangian as follows: 

C x = C + Xd^A^ + X -aX 2 + C{U X + v 2 )C (3.41) 

where C(x), acting as a Lagrange undetermined multiplier, is the new scalar variable conjugate to the ghost variable 
C(x). At present, we are ready to formulate the quantization of the a — lo model. As we learn from the Lagrange 
undetermined multiplier method, the dynamical and constrained variables as well as the Lagrange multiplier in the 
Lagrangian (3.41) can all be treated as free ones, varying arbitrarily. Therefore, we are allowed to use this kind of 
Lagrangian to construct the generating functional of Green's functions 

Z[J», J, rj, r?, I £] = ij V, <p, C, C, X)exp{i J d 4 x[£ A (x) 

+JM„ + Jip + rW + ifa + ZC + CZ]}. 1 ' ' 

Looking at the expression of the Lagrangian in Eq. (3.41), we see, the integral over X(x) is of Gaussian- type. Upon 
completing the calculation of this integral, we finally obtain 
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Z[J», J, % v , £, £] = ±J D(Ap, ip, C, C)exp{i J d±x[£ eff (x) 

where 

£eff =C- ^(d^A,) 2 + C{U X + v 2 )C (3.44) 

is the effective Lagrangian given in the general gauges. In the Landau gauge (a — > 0), the Lagrangian (3.44) just 
goes over to the one given in Eq. (3.33). As proved in Ref. [32], the quantization described in Eqs. (3.34)-(3.44) is 
equivalent to the quantization performed by the Faddeev-Popov approach [42] . At last of this section, we would like 
to emphasize that the ghost term in the £ e // does not couple to the other fields. But, we do not integrate it out in 
the generating functional. Keeping this term in the effective action and in the generating functional is helpful to later 
derivations of W-T identities. 

IV. WARD-TAKAHASHI IDENTITIES FOR GENERATING FUNCTION ALS 

A. BRST-transformation 

In this section, we show that the action and the generating functional in Eq. (3.43) are invariant with respect to 
a set of BRST transformations which include the infinitesimal gauge transformations of the nucleon, a meson and u> 
meson fields as well as the transformations for the ghost fields [29-31, 35]. The BRST transformations can be written 
in the form 



where 



Sip = <;Aip, 
Sip = <;Aip, 
SAj, = <;AA^ 
SC = <;AC, 
SC = 0, 
Sip = 



Aip = ig v Cip, 
Aip = —ig v Cip, 
AA^ = d^C, 



(4.1) 



(4.2) 



AC = l-d^A 



The above transformations for the nucleon, a meson and u> meson fields can directly be written out from Eqs. (2.5) 
and (3.39). The transformations for the ghost fields may be found from the stationary condition of the effective action 
under the BRST transformations for the nucleon, a meson and uu meson fields, 

SS eff = J d 4 xSC eff = J d 4 x{(SC - ^d u A v )(U x + v 2 )C + C(U X + v 2 )SC)} = 0. (4.3) 

This expression suggests that when the ghost fields undergo the transformations shown in Eqs. (4.1) and(4.2), the 
effective action is invariant. It is easy to prove that the integration measure in Eq. (3.43) is also invariant under the 
BRST-transformations owing to the Jacobian of the transformations being unity. 

B. W-T identity satisfied by the generating functionals for Green's functions 

When we make the BRST-transformations shown in Eqs. (4.1) to the generating functional in Eq. (3.43) and 
consider the invariance of the generating functional, the action and the integration measure under the transformations, 
we obtain an identity such that [29-31] 

±J V{A„, ip,f, <p.C, C) J d i x{J^(x)SA^ (x) + rj(x)SiP(x)+ , , 

Sip{x)r](x) +5C{x)i(x)}e lS ^f+ lE ^ = V ' ' 
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where E ■ $ with E — (J^, J, 77, 77, £,£) and $ = (A^, i/>, -0, C, C) stands for the external source terms appearing in 
Eq. (3.43). The Grassmann number <; contained in the BRST-transformations in Eq. (4.1) may be eliminated by 
performing a partial differentiation of Eq. (4.4) with respect to <;. As a result, we get a W-T identity as follows: 

i JV(A^,lf, v-C, C) J d 4 x{J»(x)AA^x) - r](x)A^(x) 

+AiP(x)r](x) + AC{x)t{x)}e lS 'fi+ lE -* = 0. 1 ' 

In order to represent the composite field functions AA^, At/> and Aip in Eq. (4.5) in terms of derivatives of the 
functional Z with respect to external sources, we may, as usual, construct a generalized generating functional by 
introducing new external sources (called BRST-sources later on) into the generating functional written in Eq. (3.43) 

Z[J", J,rj,r],^,^ui",v,v} = jfj D{A ll ,^,ip,C,C)exp{iS eff + iE •$ 

+iJd 4 x[u^AA fi + vAtP + A^v] ( ' 

where and v are the sources coupled to the functions AA^, and AW respectively. Obviously, and AA^ 
are anticommuting quantities, while, v, v, Aip and Aip are commuting ones. It is easy to verify that the BRST-source 
terms are invariant under the BRST- transformation because the functions AA^, A'tjj and Atp are nilpotent with 
respect to the BRST- transformations. Thus, we may start from the above generating functional to re-derive the 
W-T identity. The result is the same as given in Eq. (4.5) except that the external source terms are now extended to 
include the BRST-sources 

i / D{A^ V,0, <p.C, C) J d i x{J^(x)AA ll (x) - rj{x)A^{x) + A^(x) v (x) . , 

+AC(x)£(x)} exp{iS eff + iE-<I> + iJ d 4 x[u tl AA li + vAijj + AiJjv] =0. ( ' 1 

Clearly, Eq. (4.7) may be represented as [29-31] 

This is the W-T identity satisfied by the generating functional of full Green's functions. 

Apart from the identity in Eq. (4.8), there is another identity called ghost equation. The ghost equation may easily 
be derived by first making the translation transformation: C — > C + A in Eq. (4.6) where A is an arbitrary Grassmann 
variable, then differentiating Eq. (4.6) with respect to A and finally setting A = 0. The result is 

i / D(A„ C, C, 0, </>){£(*) +jp x + v 2 )C{x))} exp{iSeff +iE-$ , , 

+i J d 4 x[u^AAf, + vA^j + Aipv] = 1 ' 

which may be represented in the form [29-31] 

[i{x) + (U x + ^)-A-] Z [J„ ...,v] = 0. (4.10) 

On substituting into Eqs. (4.8) and (4.10) the relation Z = e lW where W denotes the generating functional of 
connected Green's functions, one may obtain a W-T identity and a ghost equation satisfied by the functional W such 
that 

= v ' 1 



and 



Z{x) + (O x + ^)* WIJ^ ...,«]= 0. (4.12) 
6£(x) 
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C. W-T identity obeyed by the generating functional for proper vertex functions 



The W-T identity in Eq. (4.11) and the ghost equation in Eq. (4.12) may be represented in terms of the generating 
functional T for proper (one-particle-irreducible) vertex functions. The functional T is usually defined by the following 
Legendre transformation [29-31] 



r[A», C, C, <p, ip,i>; Up, v, v] = W[J M , £, £,_J, rj, 77; u^, v, v] 
- J d 4 x[J^ +ZC + C£ + Jip + rjip + ifa] 

where A^^C, C, ip, -0 and ip are field variables defined by the following functional derivatives: 

M x ) = c '( x ) = c ^ = wt) 

From Eq. (4.13), it is not difficult to get the inverse transformations 
It is obvious that 

SW ST SW _ ST SW _ ST 
Su^ Suft ' Sv Sv ' Sv Sv ' 

Employing Eqs. (4.15) and (4.16), Eqs. (4.11) and (4.12) will be represented as 

r j4„r <5r sr \ sr _sr , sr sr , 1 a M 4 g 1 
= 

and 



SC{x 

When we define a new functional T in such a manner 



it follows that 



and 



SAp SA^ a 



sr _ sr 



(4.13) 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



(□ x + v z )C{x) = 0. (4.18) 



f = r + y d 4 x(^ M ) 2 , (4.19) 



51 + I0/*0M„ (4.20) 



(4.21) 



where $ = %l>,il>,u^,v and v. Upon inserting Eqs. (4.18)-(4.21) into Eq. (4.17) and noticing j^- = d^C, we arrive at 

The ghost equation represented through the functional T may be written as 

' M SiT-yr — ^ 2 C(x) = (4.23) 



5C(a;) ^^(x) 

In the Landau gauge, since v = and = 0, Eqs. (4.22) and (4.23) are respectively reduced to 
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. r ST ST ST5T 5T5T, n 

L SAf, 5u» <5-0 Sv dtp Sv J v ' 

and 

^-^=0. (4.25) 

These equations formally are the same as those for the massless Abelian gauge field theory [29-31]. 

From the W-T identities formulated in this section, we may derive various W-T identities obeyed by Green's 
functions and vertices, as will be illustrated soon later. 



V. W-T IDENTITY FOR w MESON PROPAGATOR AND RENORMALIZATION OF THE 

PROPAGATOR 

The W-T identity satisfied by the w meson propagator can be derived from the identities shown in Eqs. (4.8) and 
(4.10). By successive differentiations of the identity in Eq. (4.8) with respect to the sources J"{y) and £(x) and then 
setting all the sources to be zero, one may obtain 

S 2 Z S 2 Z 
9 * 8J»{x)5J»{y) l^=«=-=»=o= - a 6Z( x )8u-<(y) l^=«=-=»=o ■ ( 5 ' 1 ) 

Noticing the definitions of the lo meson and ghost particle propagators 

iD„ v (x -y)= i25 J^ 5JU{y) \j*=t=...= v =o= (0+ \T{A,(x)A„(y)}\ 0") , (5.2) 

1 A {x v) = WT^n \j~=t=-=v=o= (0 + \T{C(x)C(y)}\ 0-) (5.3) 

(Here and afterward the bold letters represent the operators) and interchanging the coordinate variables and Lorentz 
indices, Eq. (5.1) may be written as 

d^D^x -y)=ia (0+ |T*{AA„(x)C( W )}| 0") 
= -ad* A(x-y) 

where T* symbolizes the covariant time-ordering product and the definition of AA M given in Eq. (4.2) has been 
considered. Similarly, when taking the derivative of Eq. (4.10) with respect to the source £(y) and then letting all 
the sources to vanish, we get 

{a x + v 2 )A(x-y) = 5 4 (x-y). (5.5) 

This is the equation obeyed by the ghost particle propagator. Differentiating Eq. (5.4) with respect to x and utilizing 
Eq. (5.5), we find 

djfflD^ix -y) = -aU x {U x + v^-H^x - y). (5.6) 

This just is the W-T identity satisfied by the full lo meson propagator. 
By the Fourier transformation 

D^(x -y) = J d 4 xD^(k)e-^*~y\ (5.7) 

Eq. (5.6) becomes 

k »k»D» v {k) = - w -^. (5.8) 
The propagator D tlv (k) may be decomposed into a transverse part and a longitudinal part: 
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JV(fc) = D T (k 2 )(g, v - ^) + D L (k 2 )^. (5.9) 
Substitution of Eq. (5.9) into Eq. (5.8) gives rise to 

In comparison of the above expressions with the free propagator which was given in the indefinite-metric approach 
previously [29] and may easily be derived from the generating functional in Eq. (3.43) by the perturbation method 
[32] 

n(o)(]A— ( 9^ - k^K/k 2 ak^k u /k 2 

one can see that the longitudinal parts in Eqs. (5.9) and (5.11) are the same, implying that the longitudinal part of 
the oj meson propagator does not undergo renormalization. 

To derive the expression of the function Drik 2 ), it is convenient to start from the Dyson equation satisfied by the 
full u> meson propagator [29-31, 36] 

D, v (k) = D$(k) + Df x {k)U^(k)D T Ak) (5.12) 

where n AT (fc) stands for the vacuum polarization operator, or say, the self-energy operator of the u> meson. Contraction 
of Eq. (5.12) with and use of the expressions in Eqs. (5.9)-(5.11) yield [1] 

fc A n AT (fc)=0. (5.13) 

This is the W-T identity obeyed by the vacuum polarization operator which is a consequence of the gauge symmetry 
of the theory. The above identity indicates that the operator IP" (k) is transverse and therefore can be written in the 
form 

11"" (fc) = (k 2 9fW - k^K)U(k 2 ) (5.14) 

where II(fc 2 ) is a scalar function characterizing the vacuum polarization. With the above representation, it is easy to 
find from Eq. (5.12) that 

= - Pti + n^)]-^ - (5 ' 15) 

Thus, the full propagator in Eq. (5.9) can be written as 

D (k\ - f 9 ^ v - k ^ kv l k2 I ak nk u /k 2 

UpuW- i fc 2[ 1+rI (P)]- TO 2 +i€ + k 2_ u 2 +ie i- ^ lb > 

When the gauge parameter a is taken to be and 1, we obtain the propagators given in the Landau gauge and in 
the Feynman gauge respectively. When the a tends to infinity, we have the propagator given in the so-called unitary 
gauge. In the lowest order perturbative approximation, the latter propagator is of the form [1-5, 29-31] 

pW(fc) = - g %~ W< (517) 

^ w k 2 - ml + it ' 

This propagator was originally derived in the canonical quantization from the vacuum expectation value of the time- 
ordered product of the transverse field operators, iD$(x — y) = (0 |T{At ai (x)At^(j/)}| 0) and by making use of 
the Fourier representation of the transverse field operator At^x) in which the lo meson momentum k is put on the 
mass-shell, k 2 — m 2 , so that the propagator in Eq. (5.17) is transverse only for this momentum [29, 36]. However, 
due to the on-shcll property of the momentum in Eq. (5.17), when evaluating the contraction k ll D^'}{k), as we see, 
there appears an indefinite result since the numerator and the denominator in Eq. (5.17) all come to zero. Especially, 
in the zero-mass limit, there is a serious contradiction that the vector field part of the Lagrangian in Eq. (2.1) is 
converted to the massless one, but the propagator in Eq. (5.17) does not and is of an awful singularity. In contrast, 
for the propagator in Eq. (5.11), the momentum is off-shell, k 2 ^ m 2 ,. Therefore, for the transverse part of the 
propagator (or say, the propagator given in the Laudau gauge), we have k tl D^ u (k) = 0, showing a definite result. 
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Moreover, in the calculation of a loop diagram involving internal to meson lines in which the momentum of the u> 
meson line is off-shell, it is necessary to use the propagator in Eq. (5.11). In particular, the good ultraviolet property 
of the propagator allows us to perform the renormalization safely (In spite of whether the current conservation holds 
or not). In the zero-mass limit, the propagator in Eq. (5.11) and the vector field part of the Lagrangian in Eq. (3.44) 
simultaneously go over to the massless ones, exhibiting the logical consistency of the theory. 

Now let us discuss renormalization of the u> meson propagator. According to the conventional procedure of renor- 
malization, the divergence included in the functions II(fc 2 ) may be subtracted at a renormalization point, say, k 2 = ji 2 
where (i may be real or imaginary, corresponding to the subtraction point being timelike or spacelike, 

n(fc 2 ) = n( M 2 ) + n c (fc 2 ) (5.18) 

where II(/^ 2 ) and II c (fc 2 ) are respectively the divergent part and the finite part of the functions II(fc 2 ). The divergent 
part can be absorbed in the renormalization constant Z% which is defined as 

Z3- 1 = l + n(/x 2 ). (5.19) 

With this definition, on inserting Eq. (5.18) into Eq. (5.16), the uu meson propagator will be renormalized as 

£>„„(*) = Z 3 D Rflv {k) (5.20) 

where 

n (u\ — r 9» v Z k^k v /k 2 a R k^k v /k 2 1 oi\ 

VRWW- i k 2_ {mR) 2 +IlR{k 2 )+l£ + ^- V 2+ieY [ ' 

is the renormalized propagator in which is the renormalized mass, a R the renormalized gauge parameter and 
Un(k 2 ) denotes the finite correction coming from the loop diagrams. They are defined as 

™5 = \fZwi u ,a R = Z^a,IV R {k 2 ) = Z 3 k 2 U c (k 2 ). (5.22) 

It is noted that the spurious mass v is a renomalization-invariant quantity, v 2 = am 2 = ai?ra 2 K = v 2 R . Especially, at 
the renormalization point, H R (fi 2 ) = 0, as seen from Eqs. (5.18) and (5.22). In this case, we have a renornalization 
boundary condition such that 

n 1 7 \ | ( gp,v — k^ky/k 2 a^k^kjjjk 2 

D R ,Ak) -{ k 2_ [m n Y+ie + k 2_ v 2 + ie } (5-23) 

which is of the form of free propagator except that the parameters are replaced by the renormalized ones. 



VI. W-T IDENTITY FOR THE VECTORIAL VERTEX FUNCTION AND RENORMALIZATION OF 

THE VERTEX AND THE NUCLEON PROPAGATOR 

A. W-T identity for the vectorial vertex function 

The W-T identity for the vectorial vertex (nucleon-nucleon-u; meson vertex) can be derived by differentiating the 
identity in Eq. (4.8) with respect to the sources £,{x),rj(y) and rj(z) and then turning off all the sources. The result 
derived, written in the operator form, is 

5 L^(0+|T{A,(^(j,)^)}|0-) 

= i (0+ \T{C(x)i>(y)C(zyp(z)}\ 0") (6.1) 
+i(O*-\T{C{x)C(y)1>{y)1>(z)}\0-) 

where the definitions written in Eq. (4.2) have been used. Similarly, by differentiating the ghost equation in Eq. 
(4.10) with respect to the sources £(y),r](y) and r](z) and then letting the sources vanishing, one may derive the 
following equation: 

5\x-y)(0+\T{^(y)^(z)}\_0-) 

= -i(D x + v 2 ) (0+ \T{C(x)C(y)^(y)ij(z)}\ 0") . 
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Here it is noted that since there is no coupling between the ghost field and the fermion field, the two fields can not 
construct a connected Green's function. Therefore, we can write 

(o+|T{C(x)C(jMy)^)}|o-) 

= (0+ \T{C(x)C(y)}\ 0-) (0+ \T{i,{y)^{z)}\ 0") (6.3) 
= - A (x - y)S F {y - z) 

where 

(0+ \T{il)(xj$(y)}\ 0") = iS F {x - y) (6.4) 

is the nucleon propagator. It is easy to verify that once Eq. (6.3) is substituted into the right hand side of Eq. (6.2) 
and applying Eq. (5.5), we just obtain the expression on the left hand side of Eq. (6.2). Acting on both sides of Eq. 
(6.1) with the operator U x + v 2 and employing the decomposition in Eq. (6.3) and the ghost equation in Eq. (5.5), 
we obtain the following W-T identity: 

— (D x + v 2 )dHG^x,y, z) = i[S\x -y)- S\x - z)]S F (y - z) (6.5) 
ag v 

where 

G^x, y, z) = (0+ |T{A M (a#(y)^(*)}| 0") (6.6) 

is the three-point Green's function which is connected. This Green' function has the following one-particle irreducible 
decomposition [29-31] 

G^x, y,z) = J d 4 x'd 4 y'd 4 z'iD^(x - x')iS F {y - y')iT v {x', y', z')iS F (z' - z) (6.7) 

in which T u (x' ,y' , z') is the vectorial proper vertex. On inserting Eq. (6.7) into Eq. (6.5), through the Fourier 
transformation, we get in the momentum space that 



^(k 2 - v 2 )k»D^{k)S F {p)T»{ Pl q, k)S F {q) 
= {2n) i 8\k+ P -q)[S F {q)-S F {p)}. 



4,4,, , „ „„c_^ o_/„m ( 6 - 8 ) 



Considering that the energy-momentum conservation holds at the vertex, we can write 

r„(p,g,fc) = (2n)H 4 (k+p-q)ig v [ lll + A^p,q)} (6.9) 

With this representation and noticing 

(k 2 - v 2 )k»D^(k) = -ak v = a{p - q) v , (6.10) 

one may obtain from Eq. (6.8) that 

{p-qY[l li + K(P^)]=S F 1 {p)-S F 1 {q). (6.11) 

It is well-known that the general expression of the nucleon propagator S F (p) can be found from the Dyson equation 
satisfied by the propagator, as was similarly done in Eqs. (5.12)-(5.16) for the ui meson propagator. The inverse of 
the propagator can be written as 

S F 1 {p)=v-M-^{p) (6.12) 

where p =7 M p M and E(p) is the nucleon self-energy. Noticing this expression, when we differentiate the both sides of 
Eq. (6.11) with respect to p^ and then set q = p, it is found that [1] 

A,(P,P) = -^. (6-13) 

This is the W-T identity which establishes the relation between the vectorial proper vertex and the nucleon self-energy. 

It is interesting to note that the above identity determines the subtraction fashion of the nucleon self-energy. As 
one knows, the divergence in the vertex A^(p, q) may be subtracted at the renormalization point p, in such a way 
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A M (p,p) = L 7 M + A» (6.14) 

where 

L = A„(p,p) | P=M (6.15) 

is a divergent constant. Substituting Eq. (6.14) into Eq. (6.13) and then integrating Eq. (6.13) over p M from p° to 
Pft, we get 

E(p) - S( M ) - L(p- M ) - (p-m)C(p 2 ) (6.16) 

where we have chosen the p° to meet 7 M p° = /U and set the integral J p dp^A^(p 2 ) = (p— fi)C(p 2 ) with the consideration 
that when p = p°, the integral vanishes and A£(p) is finite, satisfying the boundary condition A°(p) | P=M = so that 
the C{p 2 ) is also finite, having the boundary condition C(m 2 ) = 0. When the divergent constants S(/i) and L are set 
to be 



E(/x) = A (6.17) 

and 

L = -B, (6.18) 

Eq. (6.16) will be written in the form 

Z(p) = A+(p-rilB-C( P 2 )} (6.19) 
This is the formula that gives the uniquely correct way for the subtraction of the nucleon self-energy. 



B. Renormalization of the nucleon propagator and the vectorial vertex 



Based on the representation of the self-energy in Eq. (6.19), the full nucleon propagator may be written as 



1 

p(l — B) — (M + A~ fiB) + (p-p)C{p 2 ) ' 



MP) = —i m TH , * ..m , 7Z ,.^ v (6-20) 



With the renormalization constant of the nucleon propagator defined by 

Z- 1 = 1 — B, (6.21) 

the nucleon propagator will be renormalized as 

S F ( P ) = Z 2 S§(p) (6.22) 

where 

s "<<" - fct; (6 - 23 > 

is the renormalized propagator in which Mr and T,r(p) are the renormalized mass and the finite renormalization 
correction respectively. They are separately represented in the following: 

M R = Z^M (6.24) 

where Zm is the nucleon mass renormalization constant defined by 



and 

S fl (p) = -Z 2 (p-ti)C( P 2 ) (6.26) 
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with the boundary condition | = which leads to the boundary condition of nucleon propagator like this 

1 

TP-Mr + it 



Clearly this propagator is formally the same as the free propagator. 

We would like to mention here the renormalization of the vertex function defined by V ^ (p, q) = 7^ + A M (p, g) . In 
view of the subtraction in Eq. (6.14) and the following definition of vertex renormalization constant Z\\ 

Zf 1 = 1 + L, (6.28) 

the vertex function will be rcnormalized as 

f „(p, g) = Z^f >, g) - Zf 1 ^ + A>, g)] (6.29) 

where 

A«(p,g) = Z 1 A-(p,g). (6.30) 

is the finite renormalization correction to the vertex. From the boundary condition h°{p) | P=M = mentioned before, 
it follows A^(p, g) | p=q=A j= by which we have 

f«(p,g) | P=q=M =7M- (6-31) 

This just is the boundary condition for the renormalized vertex function T^(p,q) under which the vertex is of the 
form of the bare vertex. In particular, from Eqs. (6.18), (6.21) and (6.28), it is clear to see that 

Z 2 = Zl (6.32) 

This is the Ward identity satisfied by the nucleon propagator renormalization constant and the vertex one. 

At last, it is pointed out that the identities shown in Eqs. (6.13) and (6.32) and the subtraction represented in 
Eq. (6.19) formally are the same as those in QED because they are all the consequence of U(l) gauge symmetry. 
Originally, the identities mentioned above follow from the current conservation. This result is natural because the 
current conservation, as generally proved in the Gauge Field Theory [29-31], can be derived from the global gauge 
symmetry or the local gauge symmetry. It would be emphasized that the aforementioned identities hold not only for 
the case where the u> meson is considered only, but also for the general case that the u> meson and the a meson arc 
taken into account together. We take one-loop diagrams to illustrate this point. The one-loop nucleon self-energy in 
the a — co model is represented in Fig. (la) and Fig. (lb). The one-loop vectorial vertex is shown in Fig. (2a) and 
Fig. (2b). In the figures, the solid line designates the free nucleon propagator iSp (p) represented in Eq. (6.27), the 

wavy line stands for the free u> meson propagator iD^}{k) written in Eq. (5.11), the dashed line denotes the free a 
meson propagator which is of the form 

iA(q) = - L— _ (6.33) 

qZ _ m z _|_ l£ 

the bare vectorial vertex and the bare scalar vertex r^ ) are 

r (0) - in -v 

rUZ-^ (6-34) 

The above Feynman rules are easily derived from the generating functional in Eq. (3.43) by the perturbation method. 
Applying the Feynman rules, for the one-loop nucleon self-energy defined by — iS(p), we have the following expression 

E(p) - i J d^gl^sfip - k)^D^(k) + glsfip - k) A<°) (k)} (6.35) 

where the first term and the second one are respectively given by Fig. (la) and Fig. (lb). While, for the one-loop 
vectorial vertex, in accordance with the definition in Eq. (6.9), we have 
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+fl2S^( 9 -*)7^ ) (p-*)A( )(*)]. ( j 

where the first and second terms are given by Fig. (2a) and Fig. (2b) respectively. By making use of the derivative 

-^S F ( P -k) = S F (p - k) ltl S F (p - k), (6.37) 

it is easy to find that the identity in Eq. (6.13) holds. Thus, the correctness of the W-T identity in Eq. (6.13) which 
follows from the U(l) gauge symmetry of the model is verified by the perturbative calculation. The identities in Eqs. 
(6.13) and (6.32) will be helpful to facilitate calculations of the renomalization of the a — uj model. 



VII. RENORMALIZATION OF THE a MESON PROPAGATOR AND THE SCALAR VERTEX 

For later convenience, it is necessary to give a general description for the renormalization of the a meson propagator 
and the scalar vertex. We start from the Dyson equation satisfied by the a meson full propagator i A (q) 

A(q) = A'"' (q) + A*"' (q)n(q) A (q) (7.1) 

where A* (q) is the a meson free propagator shown in Eq. (6.33) and — ifl(q) represents the a meson self-energy. 
From Eq. (7.1) it may be solved that 

A («) = ~ 2 1 -nTTX ■ ■ ( 7 - 2 ) 

q z — m z a — U{q) + is 

The self-energy can be Lorentz-covariantly decomposed into 

n(q) = Q 1 (q 2 )q 2 + n 2 (q 2 )ml. (7.3) 

The divergence in the 0(g) can be subtracted at the renormalization point \x in such a way 



Q 1 (q 2 ) = Q 1 (^) + n c 1 (q 2 ) 
2 (g 2 )-0 2 (/i 2 ) + ^(g 2 ). 



(7.4) 



On substituting Eq. (7.4) in Eq. (7.2), the propagator A(q) will be renormalized as 

A(q) = Z' 3 Aft (q) (7.5) 

where 

Z'f 1 = 1- Oi(^ 2 ) (7.6) 
is the renormalization constant of the a meson propagator, 

q 2 - {m%y - tt R (q) + is 

is the renormalized propagator in which 

= Z° m x m a (7.8) 

is the renormalized a meson mass with 

Z? n = [Zi(l + {l 2 ( P 2 ))]-i (7.9) 

being the renormalization constant of the a meson mass and 

n R (q) - Z'atfSlltf) + mln c 2 (q 2 )} (7.10) 

is the finite correction to the renormalized propagator. Obviously, the fln(q) has the boundary condition 
flftiq) \ q 2 =fi 2= which leads to the boundary condition of the propagator as follows: 
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This propagator formally is the same as the free propagator in Eq. (6.33). 

Analogous to Eq. (6.9 ) for the vectorial vertex, the scalar vertex can be written as 

T{p, q, k) = (27r) 4 <5 4 (fc + p- q)ig s f{p, q) (7.12) 

where 

T(p,q) = l + A(p,q) (7.13) 

in which A(p, q) denotes the contribution of all higher order diagrams. When the divergence in the A(j>, q) is subtracted 
at the renormalization point \i , we have 

A(p,q) = L' + A c (p,q) (7.14) 

where 

L' = A(p,p) | P=M (7.15) 

is the divergent constant and A c (p, q) is the finite correction of the A(p, q). With the above subtraction, the vertex in 
Eq. (7.13) will be renormalized as 

f(p,q) = Z'^f R {p,q) = Z'f^l + AR&q)} (7.16) 
where Z[ is the renormalization constant of the scalar vertex defined by 

Z[- 1 = l + L'. (7.17) 

and Aji(p,q) = Z'A c (p,q) is the finite correction of the Tr(p, q) with the boundary condition A^(p, q) | p=q=AJ = 
which yields the boundary condition for the renormalized vertex Tr(p, q) as follows: 

r fl (p,«)| p = q=M =i. (7.18) 

This shows that at the renormalization point, the renormalized vertex is reduced to the form of bare vertex. 

It is interesting to note that there is an identity which holds between the nucleon self-energy and the scalar vertex. 
For example, from the expression of the one-loop scalar vertex A(p, q) shown in Fig. (3a) and Fig. (3b) 

A{p,q) =ifd i k[g 2 s S<£\q - k) A<°) (k)S F (p-k) 
+9l^Sf (q - k)Sf (p - k)YD$ (k)) 

and the following derivative 

JLs F (p-k) = S F (p-k)S F (p-k), (7.20) 

^■^ <^> 

Based on this identity and the expression in Eq. (6.19), the constant defined in Eq. (7.15) can be computed by 

Z/=«| (7.22) 
dM lp=fl dM 1 ' 

This relation will be used to simplify the calculation of renormalization of the a — cj model. 



we find 
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VIII. RENORMALIZATION GROUP EQUATION AND RENORMALIZED S-MATRIX ELEMENTS 



Suppose Fr is a renormalized quantity. In the multiplicative rcnormalization, it is related to the unrenormalized 
one F in such a way 

F = Z F F R (8.1) 

where Zp is the renormalization constant of F . The Zp and Fr are all functions of the renormalization point /j, = /ioe* 
where fj,o is a fixed renormalization point corresponding the zero value of the group parameter t. Differentiating Eq. 
(8.1) with respect to /x and noticing that the F is independent of jj,, we immediately obtain a renormalization group 
equation (RGE) satisfied by the function F R [27-30] 

M~t~~ + IfFr = (8.2) 

where 7^ is the anomalous dimension defined by 

lF = ^\nZ F . (8.3) 

We first note here that because the renormalization constant is dimcnsionless, the anomalous dimension can only 
depends on the ratio a = ^j^, ~/f= 7f(#,rj' t )j where tur and gR are the renormalized mass and coupling constant 
respectively. Next, we note that Eq. (8.2) is suitable for a physical parameter (mass or coupling constant), a 
propagator, a vertex, a wave function or some other Green function. If the function Fr stands for a renormalized 
Green function, vertex or wave function, in general, it depends explicitly not only on the scale /i, but also on 
the renormalized coupling constant gR, mass tur and gauge parameter ctR which are all functions of /j,, Fr = 
Fr{p, gR(fi), m,R(p,), aR^/i); /x) where p symbolizes all the momenta. Considering that the function Fr is homogeneous 
in the momentum and mass, it may be written, under the scaling transformation of momentum p = Xpo, as follows: 

F R {p;g R ,mR,a R ;n) = X Df F R (p ; g R , j) (8.4) 

where Dp is the canonical dimension of F. Since the renormalization point is a momentum taken to subtract the 
divergence, we may set /j, — iXoA where A = e* which is taken to be the same as in p — p X. Noticing the above 
transformation, the solution to the RGE in Eq. (8.2) can be expressed as 

FR(p;gR,m R ,aR,fi a ) = X D "e^ ~ lrlFW F R (p ; g R (X), m R (X)X~ 1 , a R (X); /j, ) (8.5) 

where 3k(A), m^(A) and a R (X) arc the effective (running) coupling constant, mass and gauge parameter, respectively. 
The solution written above describes the behavior of the function Fr under the scaling of momenta. 

How to determine the function Fr(p ; ■ ■ ■ ,^ ) on the right-hand side of Eq. (8.5) when the Fr(pq, ...) stands for 
a wave function, a propagator or a vertex? This question can be unambiguously answered in the momentum space 
subtraction scheme. Noticing that the momentum po and the renormalization point /iq are fixed, but may be chosen 
arbitrarily, we can, certainly, set /Xq = Po- With this choice, by making use of the boundary condition satisfied by the 
propagator, the vertex or the wave function as denoted in Eqs. (5.23), (6.27). (6.31), (7.11) and (7.18), we may write 

FR(po;gR,m R ,aR,Li) \ P 2 = ^= F R 0) (p ;g R ,m R ,a R ) (8.6) 

where the function F R (p; g R , tur, a R ) is of the form of a free propagator, a bare vertex (if the vertex is fundamental, 
i.e., follows directly from the interaction Lagrangian) or a free wave function and independent of the rcnormalization 
point (see the examples given in the preceding sections). In light of the boundary condition in Eq. (8.6) and considering 
the homogeneity of the function Fr as mentioned in Eq. (8.4), one can write 

X Df Fr( Pq ; gR(X), m R(X)X- 1 ,aR(X), M o) |^=„ S = F R °\p; g R (X), m*(A), a R (X)) (8.7) 

where the renormalized coupling constant, mass and gauge parameter in the function F R °\p, ...) become the running 
ones. With the expression given in Eq. (8.7), Eq. (8.5) will finally be written in the form [20] 

F R (p;g R ,m R ,a R ) ^ eX ^ f{X) F ( R a) (p;g R (X),m R (X),aR(X)). (8.8) 
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For a gauge field theory, the anomalous dimensions shown in Eq. (8.8) are all cancelled out in S-matrix elements. 
To show this point more specifically, let us take the two-nucleon scattering taking place via ui meson exchanges 
as an example. The exact matrix element for the two-nucleon scattering can be written out from the well-known 
reduction formula which establishes the relation between an on-mass-shell S-matrix element and the corresponding 
off- mass-shell connected Green's function [29, 36]. A connected Green's function may conveniently be derived from the 
generating functional W for connected Green's functions as mentioned in Sec. IV. For the nucleon-nucleon scattering, 
the S-matrix element is related to the following four-point connected Green's function 

G c (x u x 2 ;y u y 2 ) = (0+ \T[i;(x 1 )ij(x 2 )^(y 1 )iJ(y 2 )]\0-) c (8.9) 

where the subscript C marks the connectivity of the Green's function. According to the familiar procedure of 
irreducible decomposition [29-31], the connected Green's function in Eq. (8.9) can be decomposed into three one- 
particle irreducible ones as represented graphically in Figs. (4a)- (4c). In each of the diagrams, there are four external 
legs which represent the full of-mass-shell nucleon propagators. These propagators will be converted to the full on- 
mass-shell nucleon wave functions by the reduction formula. The shaded blobs in the diagrams stand for the exact 
proper (one-particle irreducible) vertices. Let us first concentrate our attention on the diagrams in Figs. (4a) and 
(4b). These two diagrams represent the two-nucleon scattering taking place in the t-channcl via a ui meson exchange. 
The wavy line with a white blob in the figures denotes the full lo meson propagator. Considering the well-known 
fact that a S-matrix element expressed in terms of unrenormalized quantities is equal to that represented by the 
corresponding renormalized quantities, the scattering amplitude given by Figs. (4a) and (4b) may be written as [20, 
29, 36] 



T$ = *Zfl(9i)r&(«i,Pi)u£(Pi)^£ (k)u p R WR (92,f> 2 )<(p 2 ) 

-<(92)rS(«a,PiK(Pi)i £, u,(W(9i)r^(gi,p 2 )«i(P2) 



(8.10) 



where k = q\ — p\ = p 2 — q 2 , u R {p), T R (qi,pi) and iDj} u (k) are the nucleon wave function, the proper vectorial vertex 
and the uj meson propagator respectively which are all renormalized. The renormalization constants of the wave 
function, the propagator and the vertex are denoted by \ /f Z 2 , Z 3 and Zp respectively. The constant Zr is defined by 

Z r = Z^ZP (8.11) 

because the vertex in Eq. (8.10) is now defined by containing a vectorial coupling constant g^' multiplied with an 
imaginary number i in it. The renormalized coupling constant is defined as 

0? = ^^9v (8.12) 

On the basis of the formula given in Eq. (8.8), the renormalized nucleon wave function, meson propagator and 
vertex can be represented in the forms as shown separately in the following. For the nucleon wave function, we have 



u£(p) = e£^- (A) U <°>(p) (8.13) 



where 



is the renormalized wave function which formally is the same as the free wave function, but, the M R (X) in it is a 
running mass and 

7» = ^InZ a (8.15) 

is the anomalous dimension of the nucleon wave function. 
For the renormalized uj meson propagator, we can write 

iD^(k) = e^^ W iD^°\k) (8.16) 

where 
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is the free propagator with to^(A) and a R (X) in it being the running ui meson mass and gauge parameter and 

73(A) =//A In Z 3 (8.18) 

is the anomalous dimension of the propagator. 
For the renormalized vertex, it reads 

rS(*,Pi) = ef> ^ {X) T^( qi , Pj ) (8.19) 

where i, j = 1, 2, 

T^{ qi , Pj )=ig*{\)Y (8.20) 
is the bare vertex in which g^(X) is the running coupling constant and 

Jv(\) = fi-^lnZ v = -n^-]nZ 2 - \^\nZ z (8.21) 
an an 2 dii 

is the anomalous dimension of the vertex here the relation in Eq. (8.11) has been used. 

Upon substituting Eqs. (8.13), (8.16) and (8.19) into Eq. (8.10) and noticing Eqs. (8.15), (8.18) and (8.21), we 
find that the anomalous dimensions in the S-matrix element are all cancelled out with each other. As a result of the 
cancellation, we arrive at 

This expression clearly shows that the exact t-channel S-matrix element of the two-nucleon scattering can be repre- 
sented in the form as given by the tree diagrams shown in Figs. (5a) and (5b) provided that all the physical parameters 
in the matrix elements are replaced by their effective (running) ones. 

Next, let us turn to the diagram in Fig. (4c). In the diagram, the shaded blob with four amputated external legs 
represents the nucleon four-line proper vertex. The direct term of the scattering amplitude given by Fig. (4c) can be 
represented in terms of the renormalized quantities as follows: 

T { ff = u R (q 1 )u R (q 2 )T R (p 1 ,p 2 -,q uq2 )u R (p 1 )u R (p^ (8.23) 

In accordance with Eq. (8.8), the renormalized vertex T R ( P i,p 2 ; Qi, 92) is of the form 

r R ( Pl ,p 2 ; qi ,q 2 ) = eli ^ r{X) rf(p 1 ,p 2 ;q 1 ,q 2 ) (8.24) 

where 7r(A) is the anomalous dimension of the vertex which is determined by the renomalization constant Zr = Z^ 2 
(which is the inverse of the renormalization constant of the nucleon four-point Green's function). According to Eq. 
(8.3), 

7r (A) = A^lnZr = -2/^lnZ 2 (8.25) 

Substituting Eqs. (8.13) and (8.24) into Eq. (8.23) and noticing Eqs. (8.15) and (8.25), we also find that the 
anomalous dimensions are all cancelled out. Thus, we have 

^^^^i^fe)^^,^;?!,^)^^)^^) (8.23) 

As mentioned in Eq. (8.8), the vertex r^\pi,p 2 ; qi, q 2 ) is given at p 2 = qf = jj, 2 (i = 1,2) and the physical parameters 
in it are all running ones. Since the unrenormalized vertex T(pi,p 2 ; q\, q 2 ) is not fundamental, it has a complicated 

structure, containing a series of tree and loop diagrams [36]. The expression of the vertex T R (pi,p 2 ;qi,q 2 ) can be 
determined by the perturbation method. Unlike the loop-expansion, the perturbation series of the S-matrix usually 
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is expanded in powers of the coupling constant g v . The lowest order approximation of the vertex T(pi,p2; q\, qi) 
is of the order of g A and contains two terms which are given by the truncated subdiagrams ( the box and crossed 
box diagrams) obtained from Figs. (6e) and (6f) by amputating the external lines. The scattering amplitude given 
by the tree diagrams in Figs. (6e) and (6f) and their corresponding exchanged counterparts are convergent. We 
may dress these diagrams by replacing the free wave functions, the free propagators and the bare vertices with the 
exact wave functions, the full propagators and the rigorous proper vertices. In this way, we obtain a series of loop 
diagrams. As mentioned before, the dressed wave functions, propagators and vertices in the S-matrix clement can all 
be replaced by the renormalized ones. Therefore, they can be expressed in the forms as given in Eqs. (8.13), (8.16) 
and (8.19). Due to the cancellation of the anomalous dimensions, we will obtain an expression of the renormalized 
scattering amplitude which is formally the same as that written from the tree diagrams in Figs. (6e) and (6f) and 
their exchanged ones. For this reason, the tree diagrams are called skeletons of the dressed diagrams. There are a 
series of skeleton diagrams (or called tree diagrams) of the vertex T(j>\,p2\ qi, 92) such as the ladder diagrams and 
some others. But, in practical calculations, it is only feasible to consider the skeleton diagrams given in lower order 
perturbative approximations. We would like to stress that the skeleton diagrams can all be dressed. The S-matrix 
elements given by the dressed diagrams can be written out from the corresponding skeleton (tree) diagrams provided 
that the physical parameters are replaced by the solutions of their RGEs. For other S-matrix elements representing 
other processes, the conclusion is completely the same. It is noted that a S-matrix element evaluated by the a — lo 
model is independent of the gauge parameter as illustrated in the appendix A in the one-loop approximation (This 
is the so-called gauge-independence of S-matrix which is implied by the unitarity of S-matrix elements). This fact 
indicates that the task of renormalization for the a — lo model is reduced to find the effective coupling constants and 
the effective masses by solving their RGEs. These effective quantities completely describe the effect of higher order 
loop corrections. As an illustration, the effective coupling constants and the effective masses given in the one-loop 
approximation will be derived and discussed in detail in the next section. 



IX. ONE-LOOP EFFECTIVE COUPLING CONSTANTS AND MASSES 

For the renornalization of the a — lo model, we need to derive the effective vectorial coupling constant, the effective 
scalar coupling constant, the effective nucleon mass and the effective lo meson and a meson masses. The one-loop 
expressions of these effective quantities will be derived and discussed in the following subsections. 



A. Effective vectorial coupling constant 

The RGE for the renormalized vectorial coupling constant g^ which appears in the vectorial vertex may be imme- 
diately written out from Eq. (8.2) by setting F = g v , 

A^^)+7>).9*(M)=0. (9.1) 

In view of the definition shown in Eq. (8.1) and the relation in Eq. (6.32), the renormalization constant defined in 
Eq. (8.12) will be represented as 

Z v g = ~~r = ZP (9.2) 
Z 2 Z~i 

According to the definition in Eq. (8.3), the anomalous dimension jg(fi) in Eq. (9.1) can be calculated by 

7" = lim/x-^-ln^ = -I lim^lnZa (9.3) 

where e = 2 — ^ with n being the dimension of the space in which the regularization is performed. Based on the 
definition denoted in Eq. (5.19), the renormalization constant Z% will be given by the subtraction of the lo meson 
vacuum polarization (self-energy) operator — iH fll ,(k). From the one-loop diagram represented in Fig. (7), one may 
write [1, 29] 
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where the factor 2 comes from the fact that the fermion loop can be formed by both of the proton and neutron loops. 
By the dimensional regularization [28-30, 43], the divergent integral shown above is easily calculated. Then, from the 
definitions in Eqs. (5.14) and (5.19) , it is not difficult to find that in the n-dimensional space, the renormalization 
constant Z 3 is expressed as 

Z3 = 1 - n( <" 2) 

where m g is a mass introduced to make the coupling constant to be dimcnsionless in the n-dimcnsional space. It is 
noted here that the factors (iirm 2 ) 5 and T(l + e) may all be set to unity because they do not give an effect on the 
anomalous dimension when we set e — > in the final step of the calculation for the anomalous dimension. Inserting 
Eq. (9.5) into Eq. (9.3), it can be found that 

where a = — . In this expression, the coupling constant g v and the nuclcon mass M are unrcnormalizcd. In the 

approximation of order g 2 , they can be replaced by the renormalized ones <?^and Mr because in this approximation, 
as pointed out in the previous literature [28] , the lowest order approximation of the relation between the g v (M) and 
the g^(Mfi) is only necessary to be taken into account. Furthermore, when we introduce the scaling variable A defined 
by [i — /ioA for the renormalization point and set = Mr (This can always be done since the /iq is fixed, but may 
be chosen at will. The above choice amounts to take the renormalization scale parameter to be the nucleon mass), 
we have a = = j. Thus, with the ^ expressed in Eq. (9.6) and noticing /i^ = A^, Eq. (9.1) may be rewritten 
in the form 

} dg*(\) k fl (A)] 3 



where 



in which 



{ 



7W-t coth "^ ^ A > 2. 



Upon substituting Eqs. (9.8) and (9.9) into Eq. (9.7) and then integrating Eq. (9.7) by applying the familiar 
integration formulas, the effective (running) coupling constant will be found to be 



where a v R (X) = [g v R (X)] 2 / in , a v R = a v R (l) which is the coupling constant that should be determined by experiment 
and 

G v{ X ) = Ji x F s( A ) fan) 
= 2 + V3^-^ + (l + £)i^(A) 1 • ; 

<p(\) = VA2~41n i(A + VA^4) 

- r - V4~A 2 cos" 1 f , if A < 2, (9.12) 
~ 1 VP - Icosh" 1 1, if A > 2. 

As mentioned before, the variable A is also the scaling parameter of momenta, p — Apo, and it is convenient to 
put po 2 = ^o 2 s ° as t° apply the boundary condition. Thus, owing to the choice fi — Mr, we have po 2 = Mr 2 and 



in which 
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A = {ttt) 2 • In this case, it is apparent that when A = 1, Eq. (9.10) will be reduced to the result given on the nucleon 

mass shell with the value a v R {\) = a R . The behaviors of the effective coupling constants obtained in the timelike and 
spacelike momentum space subtractions are separately represented in Figs. (8) and (9). The figures show that the 
effective coupling constant a R (X) has a singularity Ao (the Landau pole). The position of the pole strongly depends on 
the parameter a v R . By our numerical test, we find, when the a v R is getting smaller and smaller, the Ao is getting larger 
and larger. If the a R goes to zero, the Ao approaches a value near infinity, similar to the case of QED. While, when 
the a v R is getting larger and larger, the Ao moves toward unity; but it can not arrive at unity because ce R (l) = a v R . In 
the region [0,1] of A, the cx v R (X) has no singularity to appear. In Fig. (8), there are two lines representing the a R (X) 
given in the timelike momentum subtraction: one is given by taking a v R = 1 and has a singularity at A ~ 1.1385; 
another is obtained by taking a v R = 0.5 and has a singularity at A ~ 1.3885. When A goes from Ao to zero, the a R (X) 
decreases and tends to zero, exhibiting an asymptotically free behavior as we met in QED. In Fig. (9), the two lines 
represent the a R (X) given in the spacelike momentum subtraction: one line is obtained by taking a v R — 1 and has a 
singularity at A ~ 26.4689; another is given by ot v R — 0.5 and has a singularity at Ao — 280.431. When A goes to 
zero, the a R (X) approaches the constant a v R . As one knows, the Landau poles mentioned above give a limitation of 
applicability of the one loop renormalization. That is to say, beyond the region [0, Ao], the a R (X) is meaningless even 
though in the limit: A — > oo, the a R (X) tends to zero from an opposite direction. In comparison of Fig. (8) with Fig. 
(9), it is clear to see that the range of applicability for the a R (X) given in the spacclikc momentum subtraction is 
much larger than the range for the a R (X) given in the timclikc momentum subtraction. 



B. Effective u> meson mass 



In Eq. (8.2), when we set F = m w , we have a RGE for the renormalized mass of uj meson such that 

From the definitions given in Eq. (8.1) and the first equality in Eq. (5.22), we find 

= = Zl (9.14) 



so that 



i*M = 7» = -^#^(A) (9.15) 



6ir 2 ~ 9 

where F£(X) was given in Eq. (9.8). With the above expression, Eq. (9.13) can be written as 



dm5 _ 2a R (X) p , m rfA 

Integrating the above equation, one gets 



™5m = mje* *"""< A » F .'I« (9.17) 

where = m^(l) is the observed u> meson mass. This is just the one-loop result of the effective u> meson mass. If 
we take the approximation a R (X) ~ a v R , the above expression becomes 

m£(A) ~ m^e^ G " (A) (9.18) 

where G V (X) was given in Eq. (9.11). 

To have an insight into the behavior of the effective masses in Eq. (9.17), we take the m^(A) given by taking 
a R — 1 as an example. This m^(X) is shown in Fig. (10). In the figure, the solid line represents the effective 
mass obtained in the spacelike momentum subtraction and the dashed line describes the one given in the timclikc 
momentum subtraction. Comparing Fig. (10) with Figs. ( 8) and (9), we see that the both effective masses have 
the same singularities and the same scopes of applicability as the corresponding effective coupling constants a R (X). 
Particularly, the position of the singularity strongly depends on the choice of a R as the ce R (X) does. When A tends 
to zero, the m^(A) for the spacelike momentum approaches a nonvanishing value near the m^, while, the m^(A) for 
the timelike momentum goes to zero. 
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C. Effective nucleon mass 



as 



The RGE for the renormalizcd nucleon mass, according to Eq. (8.2), can be written 

At^MflOi) + y M (p)M R (ji) = (9.19) 

where 7m(m) is the anomalous dimension of nucleon mass. In view of Eqs. (6.25) and (8.3), we have 

7m(m) = M^ln^M = -^Ml + Z 2 [± + (1 - (9.20) 

To determine the constants A, B and Z 2 in the one-loop approximation, it is necessary to compute the nucleon self- 
energy written in Eq. (6.35). By the dimensional regularization, it is not difficult to derive from Eq. (6.35) the 
following expression 

E(p) = (p-/i)S 1 (p) + E 2 (p) (9.21) 

where 

=i(p) = (4 H + (1 - «) (i^Uo Mo «^t^+ 

g-^y^^Kl - ^)b 2 + M ( P + /i)] + M( P + M )]} - 1} (9.22) 

+73^ Jo dx ee^) e 
EaW - Jo ^ 2[( Tel } ff M1 + (1 - «) 4^0 Jo 4* 

x{ (1 :e!fe 2M + eTib^WKl - + M}} + \{M „)} (9.23) 
i gl r 1 -j t (^-i)m-m 

~T(4*-) 2 JO UX ee 3 (a;) e 

6„(x) = p 2 x(x - 1) + Af 2 x + m 2 (1 - x) 

e v (x,y) = p 2 xy(xy- 1) + AT 2 xy + m 2 [(1 - x)y + a(l - y)] (9.24) 
Q s (x) = p 2 x(x - 1) + Af 2 x + to 2 (1 - x). 

From the definition ^4 = written in Eq. (6.17) and the expressions in Eqs, (9.21)-(9.23), we find 

A = A 1 +A 2 + A 3 (9.25) 

where 

-4/'*™, (9.26) 
(4tt) 2 j en v (xy 



and 



where 



and 



in which 



A 2 = (1 - «) (4 {/o cte Jo d w { (1 to„(^)e M , q 27 x 

+T7^b)*V[(l - *V)^ + My 2 ]} + \{M — fj,)} [ ■ > 



9s f 1 , - 1)A» - m 

7TT\ ( 9 - 28 ) 



(4tt) 2 J sn s ( x y 



Q v (x) = n 2 x(x - 1) + M 2 x + ml (1 - x) 

Q„ (x, y) = [i 2 xy{xy - 1) + M 2 xy + m 2 [(1 - a) j/ + a(l - y)] (9.29) 
ft s (x) = n 2 x(x - 1) + Af 2 x + m 2 (1 - x). 
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The constant B appearing in Eqs. (6.21) and (6.25), according to Eq. (6.19), ought to be computed by 

B = (p-m)" 1 ^) - A] | P=M . 
On inserting Eqs. (9.21)-(9.29) into Eq. (9.30) and employing the formula 



one may derive 



where 



Si 



1 1 



f 

Jo 



dx 



e(b — a) 



[ax + b{l - x)Y+^ 



B — B\ + Bi + i?3 



<$? [ ^S^F - nh) x(x - 1)1(1 - 1) " 2 + 2M " 11 ' 



(9.30) 



(9.31) 



(9.32) 



(9.33) 



and 



B 2 = (1 - a) j^ijl dx ft dyy{^ + ^xy[^(2 + 7xy - 9* V) 
-2Mn(2 - 3xy)] + ni^V(l ~ wftO- ~ + M^ 3 ]} - §} 



(4 



eSl s (x) E Q, s (x 



-x(x - l)[(x- l)n 2 - Mfi]}. 



(9.34) 



(9.35) 



The terms expressed in Eqs. (9.27) and (9.34) are dependent on the gauge parameter a and look more complicated. 
However, as demonstrated in the Appendix A, the S-matrix elements evaluated in the a — ui model are gauge- 
independent. Therefore, for simplicity, these terms will not be taken into account later on. This means that we limit 
ourself to working in the Feynman gauge. 

When Eqs. (9.25) and (9.32) with the expressions given in Eqs. (9.26), (9.28), (9.33) and (9.35) are substituted 
into Eq. (9.20), noticing that Z 2 — 1 should be taken in Eq. (9.20) in the approximation of order g 2 , one may find 
an explicit expression of the anomalous dimension 7m (m) through a tedious calculation 



7m(A)= 7m ) (A)+ 7 ^(A) 



(2)/ 



(9.36) 



where 7 M ^(A) is derived from the constants in Eqs. (9.26) and (9.33), while 7 M ^(A) is given by the constants in 
Eqs. (9.28) and (9 35). The expressions of the 7 M ^(A) and 7 M ^(A) are separately described in the following. For the 
anomalous dimension 7 M ^(A), we have 



7 M ) (A) = ^S t ,(A) 

7T 



(9.37) 



in which 



where the functions £"(A) are 



S„(A) = ^(A) + £;^(A)J i (A) 



i=l 



a(A) = ^(3 + A)-^ + i^(l-A), 



(9.38) 



(9.39) 



tf(A) = £[2(l-/? 2 ) 2 -3/? 2 A 2 ](l- \) 
+ £[2(1 - (3 2 f - 5/3 2 (l - (3 2 )\ 2 - /3 2 A 4 ](1 - A), 



(9.40) 



8(A) = £[/3 2 A 2 - 3(1 - /? 2 ) 2 ](1 - A) - 5^(1 - I), 



(9.41) 
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ei(A) = ^[/3 2 A 2 -(l-/3 2 ) 2 ](l-i) 
+3^[/3 2 A 4 + 2(1 - /3 2 )/? 2 A 2 - (1 - /? 2 ) 3 ](1 - A) 



and 



a(A) = Ml-/3 4 )A 2 -(l-/? 2 ) 3 ](l-I) 



^[(1-/? 4 )A 2 
-£[(1 - /3 2 ) 4 - (1 - /3 2 ) 2 (1 + 2/3 2 )A 2 + /5 4 A 4 ](1 - A) 



(9.42) 



(9.43) 



with (3 — m^/M and the functions J, (A) are given by the integrals shown below. With defining a = (1 — (i 2 )\ 2 and 
b = [3 2 \~ 2 , we can write 



JiW = Jo dx X ( x -i)+a x +b = ln 



A 2 -i-/? 2 -y g (A) 

q(X ) — Xl-l-^ + y/rtX) 



where 



q(X) = A 4 -2(l + /3 2 )A 2 + (l-/3 2 ) 2 



MX ) = f 1 dx , * +h = ln/3 + >^±£- In " 2 - 1 -" 2 -^ , 

AK ' Jo x(x-l)+ax+b 2<JtfX) \2-l-(32 +s f^' 



h (A) - / ^^(j.^+oj+i.p 



■In 



and 



= [(1 - f3 2 ) 2 - (1 + /3 2 )A 2 ] 755 ^tt - ^4w2 - 

LV M y V r ' >P 2 <lW 9 (A) 3 / 2 A 2 -l-/3 2 + ,/ g (A) 
J 4( A ) = io ^[xfx-lHax+i] 4 



?2\ A" 



In 



A^-l-/3 2 -V9(A) 



= -(A 2 + 1 - Z? 2 )-^ - (A 2 - 1 + p 2 )—^ n . - ~ = 

V f J lW K f > q(\) 3 / 2 A 2 -l-/3 2 + ^(A) 

For the anomalous dimension 7^ (A), we can write 



in which a s R = (gf) /Air and 



where the functions £f (A) are 



72 ) (A) = ^E.(A) 



X 5 (A)=£ S (A) + ^(A)J°(A) 



6(A)— | + ^ + ^-f(l-A) 



with /3 = m a /M, 



H (A) = £[2(1 - /3 2 ) 3 - 5/? 2 (l - /3 2 )A 2 - /3 2 A 4 ](1 - A) 
-£[2(l-/3 2 ) 2 -3/3 2 A 2 ](l-±) 

£!(A) - £[/3 2 A 2 - 3(1 - /? 2 ) 2 ](1 - A) + - I), 



CI(A) = ^[(l-/3o 2 ) 2 -/?o 2 A 2 ](l-i) 
+ ^[/3 2 A 4 + 2(1 - P 2 )P 2 \ 2 (1 - /? 2 ) 3 ](1 - A) 



(9.44) 



(9.45) 



(9.46) 



(9.47) 



(9.48) 



(9.49) 
(9.50) 

(9.51) 

(9.52) 
(9.53) 
(9.54) 
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and 



a(A) = Mi-/3 2 ) 3 -(i-/3 4 )A 2 ](i-i) 

&[(1 - /? 2 ) 4 - (1 - /3 2 ) 2 (1 + 2/? 2 )A 2 + /? 4 A 4 ](1 - A) 



(9.55) 



and the functions Jf(A) formally are the same as the functions Jj(A) except that the parameter /3 in the Jj(A) is now 
replaced by /3o, 



Substituting the 7m (A) as expressed in Eqs. (9.36)-(9.56) into Eq. (9.19) and solving the equation with noticing 
[idj \i — Ad/A, we obtain 



where Mr = Mr(1) is the observed nucleon mass. The coupling constants in the above have been taken to be running 
ones. The a R (X) was given in Eq. (9.10), while, the a s R (X) will be derived in the next subsection. 

It would be emphasized that in the timelike momentum space subtraction, the scaling parameter A is real so that 
the effective nucleon mass is real, while, in the spacelike momentum space subtraction, the A is imaginary so that 
the effective nucleon mass becomes complex one. In the latter case, the A in the jm (A) should be set to be iX. 
Observing the expressions in Eqs. (9.39)-(4.48) and (8.51)-(9.55), we see that in the both subtractions, the functions 
Ji (A) are always real. The real and imaginary parts of the 7m (A) are distinguished by the real and imaginary parts 
of the functions £?(A) and £f (A). In Figs. (11)-(12), we show the behaviors of the effective nucleon masses given by 
the expression in Eq. (9.57) for which the coupling constants be taken as the constants a v R and a s R for simplicity 
of computation. The effective mass Mr(A) obtained in the timelike momentum space subtraction is exhibited in 
Fig. (11). In the figure, the dashed, dotted and solid lines represent the M R {\) given by taking (a R ,a s R ) = (0,0.5), 
(0.5,0) and (0.5,0.5) respectively. From the figure, we see that in the region [0,1] of A, the Mr(A) almost keeps a 
constant Mr. Beyond this region, if only the scalar coupling is considered, the Mr (A) increases up to a maximum at 
Ao = 4.21 and then decreases to zero rather rapidly when A — > oo. While, if the vectorial coupling enters, the Mr(A) 
increases a little when A goes to a smaller Ao and then decreases much rapidly down to zero when A varies from Ao 
to oo. The effective mass Mr(A) obtained in the spacelikc momentum space subtraction is shown in Fig. (12). In 
the figure, the real and imaginary parts of the Mr(A) are displayed separately. The dashed, solid and dotted lines 
in the figure represent the real and imaginary parts of the Mr(A) which are obtained by taking {a R ,a s R ) = (0.5,0), 
(0.5,0.5) and (1, 1) respectively. The figure indicates that in the region [0,1] of A, the real part of the Mr(A) keeps 
almost a constant equal to Mr, while, the imaginary part of the Mr(A) is almost zero. When A varies from unity to 
infinity, the real part of the Mr (A) at first increases smoothly and the imaginary part of the Mr (A) decreases, then, 
both of them drastically oscillate and damp to zero. The figure also shows that the stronger the couplings (especially, 
the vectorial coupling), the larger is the frequency of the oscillation. The appearance of the oscillation implies that 
the Mr(A) is invalid to use in the region that the oscillation appears. 



When setting F = g s in Eq. (8.2), one obtains the RGE for the renormalized scalar coupling constant which is 
included in the scalar (nucleon-nucleon-cr meson) vertex 



'0 



(9.56) 



Mr(A) = M R e~ ^ dX/XlM{X) = MjRe -/ 1 AdA /A[a? ! (A)AS„(A)+ Q ? ? (A)/27rS s (A)] 



(9.57) 



D. Effective scalar coupling constant 




(9.58) 




(9.59) 



should be calculated from the renormalization constant Z s g which is represented as 




(9.60) 
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where Z 2 , Z' 3 and Z[ were defined respectively in Eqs. (6.21), (7.6) and (7.17). According to these definitions, in the 
approximation of order g 2 and in the Feynman gauge, The Z s g will be written as 

Z s g = 1 + B 1 + B 3 - L\ - L 3 + fix (9.61) 

where B\ and B% were represented in Eqs. (9.33) and (9.35) respectively, L[ and L' 3 are the parts of the constant 
L' which can conveniently be determined by the identity in Eq. (7.22). From the identity and the representations 
written in Eqs. (9.26) and (9.28), it is easy to get 

L'i = fft = £ fo dx{^ [x(x 1)/, + 2Mx] T M_} 

L's = = S fo <k{acW + - ^ Mx ^} 

The one-loop expression of the divergent constant f2i(/z 2 ) in Eq. (9.61) can be derived from the a meson one-loop 
self-energy depicted in Fig. (13). From Fig. (13), it reads 

Sl{q) = -2ig 2 s I 7^rTr[- \- —7 ^ rr] (9.63) 

w ys J (2tt) 4 L (l-q-M + ze) (q-M + ze) J v ; 

where the factor 2 also arises from nucleon doublet. By the dimensional regularization, the above integral is easily 
calculated and expressed as 

Ctf \ — 9 2 s r 1 J 24[q 2 x(x-l)+M 2 ] 

- JO aX e[q 2 x(x-l)+M^ (Q 64) 

= fM<z 2 )<z 2 + W)™ 2 - 

which gives rise to 

(4tt) 2 y e[^ 2 a;(x - 1) + M 2 ] £ 

On substituting Eq. (9.61) into Eq. (9.59), we get 

7*(A) = 7 9 s i(A) + 7 9 S 2 (A) + 7 S S 3 (A) (9.67) 

where 7gi(A), 7* 2 (A) and 7^ 3 (A) are separately defined and described in the following. 
For the 7^1 (A), we have 

7*i(A) - - Li) = ^r"(A) (9.68) 



and 



where 



r»(A) = fjg(A) + 2f,?(A)J i (A) (9.69) 



in which 



i=l 



n »fU _ 3 3 , 



rtf(A) = £[2(1 - /? 2 ) 3 - 6(1 - /3 2 ) 2 A + (1 - /3 2 )(4 

-5/? 2 )A 2 + 9/? 2 A 3 - 3/? 2 A 4 ] 
%"(A) = -^[(1 - (i 2 ) 2 - 3(1 - (i 2 )\ + (2 - /3 2 )A 2 ] 
773(A) = 2^[/3 2 A 4 - 3/3 2 A 3 - 2(1 - /3 2 ) 2 A 2 
+3(1 - /? 2 ) 2 A - (1 - /3 2 ) 3 ] 
r,l(A) = £[(2 - /3 2 )A 4 - 3(1 - /3 4 )A 3 - (1 - /3 2 ) 2 (1 
-2/? 2 )A 2 + 3(1 - /3 2 ) 3 A - (1 - /3 2 ) 4 ] 



(9.70) 
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and the functions J,(A) were given in Eqs. (9.44)-(9.48). 
For the 7| 2 (A), we can write 



where 



7 * 2 (A) = ^(B 3 L' 3 ) = (A) (9.71) 



r?(A)= % *(A) + X> s (A)J°(A) (9.72) 



in which 



3 , l-/3 2 
"2 + 



n{{\) = £[2(1 - /? 2 ) 3 - (1 - /3 2 )(2 + 5/? 2 )A 2 ] 

,?f(A) = £[A 2 -(l-/3 2 ) 2 ] (9.73) 

^I(A) = ^[/?o 2 A 4 + (1 - $)(!_+ 2$)A 2 - (1 - /? 2 ) 31 



VtW = -Mi 1 + $W 2(1 - /3 2 ) 2 (1 + /? 2 )A 2 + (1 - /3 2 ) 4 ] 

and the functions J® (A) were defined in Eq. (9.56). 

For the 7j 3 (A), By virtue of the expression given in Eq. (9.65), one can get 



V " 8tt 2 



^(A^-fi^-^-I^A) (9.74) 



where 



Based on the anomalous dimension 7* (A) given in Eqs. (9.67), (9.68), (9.71) and (9.74), the RGE in Eq. (9.58) 
may be represented in the form 



+ P(A). g f(A)+Q(A)( 5 f(A)) 3 = (9.76) 



rfA 

where 



P(A) = ^m T v {x)t Q(A) = j_ [r? (A) _ r « (A)] (9J7) 

In the above equation, the o^(A) is a known quantity as given in Eq. (9.10). So, Eq. (9.76) is the equation used to 
determine the unknown quantity gf (A) only. To solve the nonlinear equation, we may set 

gf{X)=u(X)-^ (9.78) 

which leads Eq. (9.76) to a linear equation obeyed by the function u{\) 

- 2P(A)u(A) - 2Q(A) = (9.79) 
When setting Q{\) = 0, we obtain a homogeneous equation whose solution is 

u(A) = u(l)e 2 /i A<iAP(A) (9.80) 
In order to seek the solution of Eq. (9.79), we assume 

u(A) = v(X)e 2 Ii dXP(X) (9.81) 



32 



where v(X) is an unknown function needs to be determined from Eq. (9.79). Inserting Eq. (9.81) into Eq. (9.79), we 
get 



dv(\) 
dX 

Integrating the above equation, one obtains 



2Q(A)e~ 2 /i AdAP(A) (9.82) 



v(X) = u(l) + 2 J dXQ(X)e~ 2 Si dAP(A) (9.83) 

Combining the expressions in Eqs. (9.78), (9.81) and (9.83), the solution of Eq. (9.76) is finally given in the form 

1 + a>* R /irLr s (X) 

where a s R (X) = {g^(X)) 2 /4*7r, a s R — a s R {\) which is a parameter needed to be determined by fitting the experimental 
data, 

^(A) = e ^/ 1 AdA/Aa « (A)r " (A) (9.85) 

and 

G S (X)=J^ dX/XT s (X)K(X) (9.86) 

in which 



r s (A)=n(A)-r^(A) (9.87) 

The behaviors of the effective coupling constant a s R (X) obtained in the timelike and spacelike momentum subtrac- 
tions are separately displayed in Figs. (14) and (15). For timelike momenta, the o:fj(A) is real. In Fig. (14), there are 
four lines representing this a R (X) which are given by four groups of the parameters (a R , a R ) = (0, 1), (0.5, 1), (0, 0.2) 
and (0.5,0.2) respectively. The figure indicates that the a s R (X) has a Landau pole Ao-The poles for the four lines are 
respectively located about at Ao ~ 1.07523 (for the first lines), 1.8237 and 2.3967. Clearly, the Ao is larger if the both 
parameters (ct R ,a s R ) are smaller. In particular, the pole moves to the point near infinity when the both parameters 
(a R ,a R ) tend to zero. When A goes from Ao to zero and infinity, each a s R (X) in Fig. (14) abruptly falls to zero from 
the opposite directions. For spacelike momenta, The behavior of the effective coupling constant is described by the 
lines in Fig. (15). The two lines in Figs (15a) and (15b) are given by considering the scalar coupling only with taking 
a s R = 1 and 0.2 respectively. In this case, the ct s R (X) is real and has a pole Ao- the poles for the aforementioned 
lines are located respectively at A = 5.725 and 234.6. However, when the vectorial coupling is included, the a s R (X) 
becomes complex. In this case, the pole disappears; instead, there is a maximum to appear as shown in Figs (15c) 
and (15d). The lines in Figs (15c) and (15d) are given by taking (a R ,a s R ) = (0.5,1) and (0.5,0.2) respectively. In 
the figures, the real part of the a s R (X) is represented by the solid line and the imaginary part by the dashed line. 
From Fig. (15c), we see that the real and imaginary parts of the a s R (X) have sharp peaks corresponding to the pole 
of the upper line. The peaks are manifested more clearly by the right amplified lines. Fig. (15d) exhibits that cither 
the real part or the imaginary part varies rather smoothly due to the decrement of the parameters (a R , a s R ) and the 
larger effect of the vectorial coupling. It is noted that when A — > 0, the ct s R (X) reaches a constant, while, in the limit 
of A — > oo, the afj(A) goes to zero. In particular, The behaviors of the a R (X) tell us that the smaller the parameters 
(a R , a s R ), the larger will be the range of applicability and the a s R (X) for the spacelike momenta have a larger range of 
applicability than that for the timelike momenta. 



E. Effective a meson mass 



In accordance with Eq. (8.2), the RGE for the renormalizcd a meson mass is 

d 



A^m^(A)+ 7 -(A)m«(A)=0 (9.88) 
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where 



l° m {\)=^\nZ° m . (9.89) 

From the definitions given in Eqs. (7.6) and (7.9), it is found that in the approximation of order g 2 , the renormalization 
constant can be written as 

Z-=l-i[0 1 ( M 2 ) + 2 ( M 2 )] (9.90) 

The one-loop expressions of the divergent constants fli(fj, 2 ) and 2 (/z 2 )] were given in Eqs. (9.65) and (9.66). With 
these expressions, the renormalization constants in Eq. (9.90) can explicitly be written out. Use of this renor- 
malization constant in Eq. (9.89) yields the a meson mass anomalous dimension as follows: 

7™ = -fU(A) (9.91) 

where 

6 1 11 2 1 

n [ 6~pl + A 2 " ~ A ( A~ 2 (31 

in which Po = ^ and 



GsW = -b - -m + T2 - T(TI - 752>^( A )] ( 9 - 92 ) 



Substituting the above anomalous dimension into Eq. (9.88) and solving the equation, we obtain an expression of the 
effective a meson mass such that 

m*(A) = m?Ji *<*< A > G 'W (9 Q4) 

where — m^(l) is a mass parameter which needs to be determined by experiment. If the coupling constant a s R (X) 
is approximately taken to be a constant a s R , the integral over A can easily be calculated. In this case, we have 



TO «(A) = m^e s "^ (9.95) 



where 



2a a n 1 _i_ V3f-i 2 



S„{\) = ^[1-^ + ^(1 



(9.96) 



It is seen from Eq. (9.94) that the behavior of the effective mass to^(A) is intimately related to property of the 
effective coupling constant a s R (X). To give a view of the behavior of the to^(A), we take the m^(A) evaluated from 
Eq. (9.94) by taking {a R ,a s R ) = (1, 1) as an example. This m^(A) is shown in Fig. (16). In the figure, the dashed 
line represents the m^(X) given in the timelike momentum subtraction. This m^(A) is real and has a singularity at 
A w 1.09758 which implies that the range of applicability of the m^{\) is less than A = 1, the nucleon mass scale. 
When A tends to zero, the m^(A) approaches a value which does not deviate from the constant so much. The 
to^(A) given in the spacelike momentum subtraction is complex .The solid line in Fig. (16) represents the real part 
of the m^(A) which has a maximum near A = 5 which indicates that the m^(A) is applicable in a wide region of 
[0, 5]. Similar to the coupling constant a s R (X), when the parameters (a R , a s R ) are taken to be smaller, either the pole 
or the maximum will be shifted to the point of a large Aq. 



X. SUMMARY AND DISCUSSIONS 



In this paper, it has been argued that the a — lo model, as a constrained system, is really of U(l) local gauge 
symmetry. This enables us to quantize the a — lu model by means of the method used for quantizing the gauge field 
theory. In particular, the gauge symmetry allows us, in a consistent way, to derive various W-T identities which 



34 



provide a faithful basis for performing the renormalization of the model. As shown in Sec. V, the W-T identity 
in Eq. (5.6) satisfied by the u> meson propagator and the W-T identity in Eq. (5.13) for the vacuum polarization 
operator determine not only the structures of the propagator and the vacuum polarization operator, but also the 
renormalization fashion of the propagator and the co meson mass as shown in Eqs. (5.14), (5.16), (5.21) and (5.22). 
Especially, the W-T identity in Eq. (6.13) obeyed by the vertex gives rise to the correct manner of subtraction of the 
nuclcon self-energy as denoted in Eq. (6.19). As shown in Sec. VI, the subtraction in Eq. (6.19) leads to the correct 
representations for the renormalization constants of nuclcon propagator and nucleon mass as shown in Eqs. (6.21) 
and (6.25). Moreover, the identity in Eq. (6.13) directly yields the important relation between the renormalization 
constants Z\ and Zi as written in Eq. (6.32). This relation together with the relation in Eq. (7.22) which follows 
from the identity in Eq. (7.21) greatly simplify the calculation of the renormalization. It would be mentioned here 
that in some previous works [7, 15, 17], the subtraction based on the expression S(p) = Ap+BM was ever used. 
This subtraction gives the nucleon propagator renormalization constant as Z2 = [1 — A{ji 2 )]~ l and nucleon mass 
renormalization constant as Zm = \%i(X + -^(m 2 )] 1 which arc different from the renormalization constants written in 
Eqs. (6.21) and (6.25) and therefore the relation in Eq. (6.32) could not be fulfilled in this case. The renormalization 
of the model under consideration is performed in the mass-dependent momentum space subtraction scheme by the 
renormalization group approach. The prominent advantage of the subtraction scheme is that it naturally leads to the 
boundary conditions for the renormalized propagators, the vertices and the wave functions. The boundary conditions 
allow us to give an unique determination of the solutions to the renormalization group equations for the renormalized 
propagators, vertices and wave functions without any ambiguity. As claimed in the Introduction, we limit ourself in 
this paper to examine the renormalization of the model at zero temperature by means of the renormalization group 
method. Since the perturbative series expanded in the powers of coupling constants is chosen to be the starting point of 
this renormalization, the results of the renormalization would be different from those obtained in the study of nuclear 
matter by using the loop expansion and the spectral function methods. Hopefully, the renormalization procedure 
described in this paper will be helpful for applying the renormalization group approach to study the nuclear matter 
at finite temperature and finite density. 

The procedure of renormalization group method was demonstrated by the one-loop renormalization in this paper. 
Since the renormalization exactly respects the W-T identities, the results obtained arc faithful. Especially, the one-loop 
effective coupling constants and masses are given in the rigorous forms as they are derived from the mass-dependent 
momentum space subtraction. The subtraction scheme used is, in principle, suitable not only for high energy, but 
also for low energy, unlike the minimal subtraction scheme [28-31] which is only appropriate in the large momentum 
limit. In addition, the expressions of the one-loop effective physical quantities derived in this paper are applicable for 
the both of timelike momentum subtraction and spacelike momentum subtraction. As seen from Figs. (8)-(12) and 
(14)- (16), the behaviors of the effective quantities given in the timelike subtraction and the spacelike subtraction are 
much different from one another. In which case we should use the results given in the timelike momentum subtraction 
or in the spacelike momentum subtraction? The answer to this question depends on what process is discussed. For 
example, when we study the nucleon-nucleon scattering taking place in the t-channcl, as mentioned in the Appendix 
B, the transfer momentum in the boson propagator is spacelike. In this case, it is suitable to take the effective coupling 
constants and boson masses given in the spacelike momentum subtraction. If wc investigate the nucleon-antinucleon 
annihilation process which takes place in the s-channel, since the transfer momentum is timelike, the effective coupling 
constants and boson masses given in the timelike momentum subtraction should be used. The effect of the one-loop 
renormalization is examined by the nucleon-nucleon elastic scattering whose differential cross section given in the 
order of g 2 is described in the Appendix B and plotted in Fig. 17. In the figure, We only take the differential cross 
sections given at the laboratory kinetic energies 2} a b = 491.9-MeV and 575.5MeV r as an example. The figure shows 
that consideration of the one-loop renormalization requires the coupling constants to be smaller in order to fit the 
experimental data. This actually is a general feature of considering the renormalization effect. 

As exhibited in Sec. IX, the one-loop effective physical parameters given in the timelike momentum subtraction and 
the ones given in the spacelike momentum subtraction not only behave differently, but also have different ranges of 
applicability because the singularities of the effective quantities given by the two subtractions appear at the different 
momenta. Especially, the positions of the singularities are strongly dependent on the coupling constants a 1 ^ and 
af. The smaller the coupling constants, the larger are the ranges of applicability. It would be mentioned that 
since the propagators written in Eqs. (5.16), (6.12) and (7.2) are solved from the Dyson equations, the one-loop 
renormalization actually contains the contribution given by partially summing up a set of chain loop diagrams. Just 
due to the partial summation, as mentioned before, the coupling constants must be set to be smaller for fitting 
the experimental data of the nucleon scattering. To this end, it is natural to ask if and how the behaviors of the 
one-loop effective physical parameters can be modified by considering higher order loop renormalizations? In other 
words, if the coupling constants would be more smaller and the ranges of applicability of the renormalization would 
be enlarged when the contributions arising from more higher order loop diagrams are summed up? Obviously, this 
is an interesting problem worthy to pursuing further. In addition, we would like to address that the a — v model 
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should be viewed as a restrictive model in which the a field is introduced as a phenomenological field. Aside from 
the o — to model, there are some other models in QHD which are of a certain gauge symmetry. Especially, the model 
proposed by Sakurai in the early time [44], in our opinion, is most promising to describe the nuclear force because in 
this model, exchanges of the light mesons, pion and pho meson, dominate the strong interaction between nucleons. In 
view of the argument given in Refs. [32-34, 45], the Sakurai's model is a SU(2) gauge field theory which is not only 
gauge-invariant, but also renormalizable. Certainly, the renormalization of this model may be investigated along the 
same line as described in this paper. We will discuss this subject in the future. But, it can not be expected that a 
perturbative investigation could give an ultimate solution to the strong interaction. Just as said in Ref. [46, 47], to 
resolve the strong interaction, it is adequate to perform a nonperturbative study of the interaction kernel appearing 
in the relativistic equation whose closed expression can be derived by the procedure as described in Refs. [46, 47]. 
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XII. APPENDIX A: GAUGE-INDEPENDENCE OF S-MATRIX ELEMENTS 



The gauge- independence of S-matrix elements computed by a gauge field theory is a well-known fact. For the a — uo 
model, as argued in this paper, it actually is a U(l) gauge field theory. So, the same conclusion should hold for the 
<j — to model. To convince oneself of this fact, we take the nucleon-nucleon scattering amplitudes up to the one-loop 
approximation as examples to show that the matrix elements given by the a — oo model are surely independent of the 
gauge parameter a. The typical Feynman diagrams representing the scattering amplitudes are depicted in Figs. (5) 
and (6). In Fig. (6), only the diagrams with the internal to meson line are necessary to be considered. 

For the tree diagram in Fig. (5a), the gauge-independence of its S-matrix element is well-known. In fact, the 
S-matrix element 

Si ^u s >(q 2 )ig v ^^u r >{q 1 )iD^(k)u s {p2)igvluU r {pi) (A.l) 

where u s (p) is the free nucleon wave function can be divided into two parts according to the decomposition of free to 
meson propagator 

D^(k) = D^(k) + £>£"(*) ( A - 2 ) 

where 

D p(k) = - (A.3) 
k z — m~ + is 

which is the propagator given in the Feynman gauge and 

£>£"(*) = (! - a)D a {k 2 )W (A.4) 
which is the a— dependent part of the propagator in which 

D ( k ^ = 712 2 T ■ 1,2 T-—,- (A-5) 

(Ar - m~ + ie){k A - v 2 - + te) 

For the a— dependent part of Si, applying the energy-momentum conservation k = q\ — q 2 = P2 — Pi and Dirac 
equation (p— M)u s (p) = 0, it is found 

Sf = -(1 - a)iglu s f(q 2 )\^u rl {q 1 )u s {p2)'ku r (pi)D a {k 2 ) 

= -(! - Oi)iglu s r(q 2 ){ci 1 - q 2 )u r /(gi)u s (p 2 )(P2 - pi)u r (pi)A,(fc 2 ) (A.6) 
= 0. 

This shows that the u> meson propagator given in the Feynman gauge is sufficient to use for evaluating the tree 
diagram matrix element. In the same way, one may prove that the S-matrix element given by the tree-diagram in 
Fig. (5b) is independent of the gauge parameter as well. 
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Let us focus on the one-loop diagrams in Fig. (6) where only the direct diagrams are plotted and necessarily to 
be examined for our purpose. The gauge- independence of the matrix element of Fig. (6a) which contains a ui meson 
self-energy in it can easily be proved by using Eq. (A. 6). So, we only need to examine the gauge-independence of 
Figs. (6b)-(6f). The matrix element of Fig. (6b) with a vertex correction in it can be written as 

S 1 , = M^q h q 2 )A^ Pl ,p 2 ) (A.7) 

where 

M^(qi t q 2 ) =u s >{q 2 )ig v l v u rl (q 1 )iD u ^{k) (A. 8) 
which is gauge-independent as shown in (A. 6) and 

/d A k 
-^-^u s (p 2 )ig v -f p iS F (p 2 - k)ig v ^iS F (pi - k)ig v -f a u r (pi)iD pa (k). (A.9) 

Replacing D pa (k) by the D P t a (k) shown in (A. 4), we have the following gauge-dependent part of A fJ, (pi,p 2 ) 

A£(pi,P2) = -(1 - a)glj ^u s (p 2 )kS F (p 2 - k)^S F ( Pl - k)ku r { Pl )D a {k 2 ) (A.10) 

where k = j^k^ can be written in the form 

k=S F 1 (p l )-S F 1 (p l -k) (A.ll) 
here i = 1,2. Using this relation and Dirac equation, Eq. (A. 10) becomes 

A£(pi,P2) = -(1 - a)g 3 v u s (p 2 )^u r ( Pl )J (A.12) 

where 

J = [ -^jD a (k 2 ) = lim -t- [ — dx 

J (2tt) 4 ay ' ^o(47r) 2 i e{[a+{l-a)x]mlY y ' 

here the last equality is given by the dimensional regularization. This integral gives a divergent constant without 
containing any finite number in it. Therefore, it may completely be cancelled out by a counterterm in a renormalization 
program and gives no contribution to the renormalized S-matrix element. On the other hand, since the integral is 
independent of momentum, it would not contribute to the anomalous dimension and hence to any physical quantity. 
For the matrix element of Fig. (6c) which contains a nucleon self-energy, it may be written as 

Sl = M li (q h q 2 )B^ Pl ,p 2 ) (A.14) 

where 

/d 4 k 

j^-^u s { P2 )ig v ~f p iS F {p 2 - k)ig v ^ a iS F {p 1 + q)ig v j tl u r (pi)iD f " T (k) (A.15) 

in which the a— dependent part is of the form 

B£(Pi.P2) = -(1 - a)glj ^Lu s ( P2 )kS F ( P2 - k)kS F ( Pl + q)^u r ( Pl )D a (k 2 ) (A.16) 

By employing the relation in Eq. (A.ll) and Dirac equation, one may find 

B£(pi,P2) = (1 - a)g 3 v u s (P2hvS F ( P2 )^u r ( Pl )JZ = (A.17) 
This is because the integral in it vanishes 

J'2= J ^- 4 k»D a (k 2 ) = (A.18) 
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due to that the integrand is an odd function. Similarly, the a— dependent part of Fig. (6d) can also be proved to give 
no contribution to the S-matrix element. 

Let us turn to Figs.(6e) and (6f). The matrix elements of the both figures can be respectively represented as 

S $ = 9iJ$j*u s ,(q 2 )rS F (qi - fc)7"«r'(«i) (A 19) 
■xu s (p 2 )7 p Sf(pi + k)-fur{pi)Dw(q - k)D ua {k) 

and 

s$ = gtJ(0*uAq2)rSF(qi - k)Yu r ,{ qi ) (A 2Q) 

xu s (p2h p S F (p2 - fc)7 CT u r (pi)£) iU(7 (g - k)D up {k) 

where q = p 2 — pi = qi — q 2 . Their a— dependent parts are denoted by 5| Q and S 2a . By making use of the relation 
in (A. 11) and the relation q — k = (qi — k) — q 2 = P2 — (k + pi) as well as Dirac equation, it is easy to find 

Sl a = -2(1 - a)^(9a)7"«r'(«l)« a (p2)7M«r(Pl) / (2^ 

-(1 - afgtu sl {q 2 )^u r ,{q l )u s {p 2 )Yu r {p 1 ) J ^k^D a [(q - k) 2 ]D a (k 2 ) ( A - 21 ) 

_ _ C4 
— D 2a 

which gives Sfa + S 2a = so that the sum of S 2 and S 2 is independent of the gauge parameter. 

XIII. APPENDIX B: CROSS SECTION OF NUCLEON-NUCLEON SCATTERING 

To illustrate the effect of the renormalization described in this paper, we evaluate the cross section of the nucleon- 
nucleon elastic scattering. Here we limit ourself to first consider the cross section given in the approximation of order 
g 2 . In this approximation, only the tree diagrams denoted in Fig. (5) are concerned. From these diagrams, in the 
center of mass frame, the differential cross section is easily calculated and represented as follows: 

= i[o%T v + a 2 s T s - a v a s T vs ) (B.l) 



where S = A{p 2 + M 2 ) with p being the nucleon momentum is the squared total energy of the system, T v is contributed 
from the oj meson exchange interaction, T s is given by the a meson exchange interaction and T vs is the crossed term 
related to both of the u meson and a meson exchanges. They are separately represented as follows: 

r ~ ' 2 +(-D 1+/ ^ (B.2) 



where 



and 



(A?) 2 (A2) 2 v 1 A?A5 

R\ = M 4 + 2p 2 M 2 cos 6 + 2M 2 p 2 + 2p i + 2p A cos 4 6/2, (B.3) 

R v 2 = M 4 - 2p 2 M 2 cos 6 + 2M 2 p 2 + 2p 4 + 2p 4 sin 4 6/2, (B.4) 

Rl = 16(p 4 - M 4 ), (B.5) 

A^ = Ap 2 sin 2 0/2 + m„, (B.6) 

A^ = Ap 2 cos 2 6/2 + ml. (B.7) 

JDS JDS JDS 

7 - ' ! --fe": 1 7 T%- (B.8) 



(Af) 2 (A|)2 ^ ' Af A| ' 
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where 

Rl = 4(4p 2 sin 2 fl/2 + M 2 ), (B.9) 



and 



R s 2 = 4(4p 2 cos 2 9/2 + M 2 ) 2 , 
R s 3 = 2[2M 2 (p 2 + M 2 ) + / (sin 4 6/2 + cos 4 6/2], 
A{ = 4p 2 sin 2 6/2 + m 2 , 

A* = 4p 2 cos 2 6/2 + ml. 

jdvs I3 VS T3 VS T? v 

T - 1 I 2 ( l) 1 If 1) J 4 

" s ADAS 1 ADAS 1 V / ADAS 1 V / A V , 



AJAf A«A* v ' A- A* v ' Aj Af 
where 

i?f = 4M 2 [(Af 2 + 2p 2 ) 2 + (Af 2 + 2p 2 cos 2 <9/2) 2 ], 
i?^ = 4M 2 [(M 2 + 2p 2 ) 2 + (M 2 + 2p 2 sin 2 6/2) 2 }, 
R V 3 S = 4[2(M 2 + p 2 cos 2 6/2) 2 - M 2 (M 2 + p 2 + p 2 sin 2 6/2)] 

and 

R V 4 S = 4[2(M 2 + p 2 sin 2 6/2) 2 - M 2 (M 2 + p 2 + p 2 cos 2 6/2)]. 



B.10) 
B.ll) 
B.12) 

B.13) 
B.14) 

B.15) 
B.16) 
B.17) 

B.18) 



In the above, (6,ip) are the scattering angles, I is the isospin of the two-nucleon system, AJJ and A| (k = 1,2) are 
respectively given by the u> meson and a meson propagators and Rf (a = v,s,i = 1,2,3) are the functions coming 
from the nucleon spinor matrix elements. In Eqs. (B.2), (B.8) and (B.14), the isospin-related terms arise from the 
exchanged diagram, while, the remaining terms represent the contribution of the direct diagram. 

To consider the renormalization effect on the two-nucleon scattering, as mentioned in Sec. VIII, we may directly 
replace the coupling constant a v and a s in Eq. (B.l) by their effective ones a v R (\) and a s R (X) and the particle masses 
M, m w and m a appearing in the propagators and the functions Rf by their effective counterparts Mr(\), m^(A) and 
m^(A). To this end, it should be noted that the nucleon spinors in the S-matrix element under consideration are on 
the mass-shell, satisfying the free nucleon Dirac equation. The momenta pi in the spinors arc timclikc because they 
meet the relation p 2 = M 2 where M is real. For the renormalized spinor wave function shown in Eq. (8.14), as easily 
verified, it also satisfies the Dirac equation and the momentum in the spinor fulfills the relation p(\) 2 = M^(A) 2 
where p(X) — (E(X),p) with E(\) = (p 2 + Mr(X) 2 ) 1 ' 2 . Therefore, for the nucleon-nucleon scattering, it is adequate 
to take the effective nucleon mass given in the timelike momentum space subtraction. While, the momenta in the oj 
meson and a meson propagators, as one knows, are off-shell and spacelike in the t-channel scattering. Therefore, it 
is appropriate to take the effective coupling constants a R (X) and a s R (X) and the effective meson masses m^(A) and 
m^(A) given in the spacelike momentum space subtraction. In this section, we only examine the effect of the one- loop 
renormalization on the two-proton scattering by using the effective coupling constants and masses presented in Sec. 
IX. The differential cross sections given at the laboratory kinetic energies Ti a b — 491.9MeV and 575.5Mey are shown 
in Fig. (17a) and (17b). In the figures, the solid lines and the dashed lines represent respectively the calculated results 
with and without considering the renormalization effect and the experimental data are taken from Ref. [48] . As shown 
in Fig. (17a) and (17b), in the case without considering the renormalization effect, the theoretical parameters are 
taken to be M R = 938MeV, a% = 1.1, a£ = 1.4, raff = 782MeV and = 580MeV; while in the case of considering 
the renormalization effect, the parameters must be taken to be Mr = 938MeV, = 0.55, = 0.62, ra^ = 782MeV 
and m R = 670MeV. It is clearly seen from the figures that consideration of the renormalization has an effect that to 
fit the experimental data, the coupling constants must be set to be smaller. 
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XV. FIGURE CAPTIONS 

Fig.(l): The one-loop nucleon self-energy in the a — ui model. The solid, wavy line and dashed lines represent the 
free nucleon, ui meson and a meson propagators respectively. 

Fig. (2): The one-loop vectorial vertices in the a — ui model. The lines represent the same as in Fig. (1). 

Fig. (3): The one-loop scalar vertices in the a — u> model. The lines mark the same as in Fig. (1). 

Fig. (4): The diagrams represent the nucleon four-point one-particle irreducible Green's function. The solid line 
with a white blob represents the full nucleon propagator. The wavy line with a white blob denotes the full u meson 
propagator, The shaded blobs represent the proper vertices. 

Fig. (5): The tree diagrams of nucleon- nucleon scattering. The first two diagrams represent the interaction 
generated by the u> meson exchange. The remaining two diagrams represent the interaction mediated by the a meson 
exchange. 

Fig. (6): Some two-nucleon one-loop Feynman diagrams which are chosen to demonstrate the gauge- independence 
of the nucleon scattering matrix elements. 

Fig. (7): the one- loop diagram of the effective u meson self-energy. The solid line marks the free nucleon propagator 
and the wavy line denotes the free ui meson propagator. 

Fig. (8): The effective one-loop vectorial coupling constants a R (X) given by the timelike momentum space subtrac- 
tion. The solid line represents the coupling constant given by taking a v R = 0.5. The dashed line denotes the coupling 
constant given by a v R = 1. 

Fig. (9): The effective one-loop vectorial coupling constants a^(A) given by the spacclike momentum space subtrac- 
tion. The solid line represents the coupling constant given by taking a v R = 0.5. The dashed line denotes the coupling 
constant given by a v R = 1. 

Fig. (10): The effective one-loop uo meson masses m^(A) given by taking a v R = 1. The solid line and the dashed 
line represent the effective masses obtained in the spacelike momentum subtraction and the timelikc momentum 
respectively. 

Fig. (11): The effective one-loop nucleon masses M R (X) obtained in the timelike momentum subtraction. The solid, 
dashed and dotted lines represent the effective masses given by (a R ,a s R ) = (0.5,0), (0,0.5) and (0.5,0.5) respectively. 

Fig. (12): The effective one-loop nucleon masses M R (X) obtained in the spacelike momentum subtraction. The 
dashed, solid and dotted lines represent the effective masses given by taking {a Rl a R ) — (0.5,0), (0.5,0.5) and (1, 1) 
respectively. The upper figure describes the real part of the M R {\). Another figure shows the imaginary part of the 
Mfl(A). 

Fig. (13): The a meson one-loop self-energy. The solid line represents the free nucleon propagator and the dashed 
line denotes the free a meson propagator. 

Fig. (14): The effective one-loop scalar coupling constants a s R (\) obtained in the timelike momentum subtraction. 
The dashed and solid lines on the left represents the effective coupling constants given by (a R , a s R ) = (0, 1) and 
(0.5, 1). The dashed and solid lines on the right denote the effective coupling constants given by (a R ,a s R ) = (0,0.2) 
and (0.5,0.2). 

Fig. (15): The effective one-loop scalar coupling constants ct s R (X) obtained in the spacelike momentum subtraction. 
The four lines represent the effective coupling constants given by taking (a R ,a s R ) = (0,1), (0,0.2), (0.5,1) and 
(0.5,0.2) respectively. The solid and dashed lines denote the real parts and the imaginary parts of the coupling 
constants respectively. 

Fig. (16): The effective one-loop a meson masses m^(A) obtained by taking (a v R ,a s R ) = (1,1). The solid line 
represents the effective mass given in the spacelike momentum subtraction, The dashed line shows the real part of 
the m^(A) given in the timelike momentum subtraction. 

Fig. (17): The two proton elastic differential cross sections given at the laboratory kinetic energies Ti a b = 491.9MeV 
and 575.5MeV. The black squares show the experimental data. The solid lines represent the theoretical values 
calculated by considering the one- loop renormalization effect. The dashed lines represent the theoretical values 
without considering the one- loop renormalization effect. 
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